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ABSTRACT 

We report the first results of a study of variable point sources identified using 
multi-color time-series photometry from Sloan Digital Sky Survey (SDSS) Stripe 
82, including data from the SDSS-II Supernova Survey, over a span of nearly ten 
years (1998-2007). We construct a light-curve catalog of 221,842 point sources 
in the RA to 4 h half of Stripe 82, limited to r = 22.0 mag, that have at least 
10 detections in the ugriz bands and color errors < 0.2 mag. These sources are 
then classified by color and by cross-matching them to existing SDSS catalogs of 
interesting objects. Inhomogenous ensemble differential photometry techniques 
are used to greatly improve our sensitivity to variability and reduce contamina- 
tion by sources that appear variable due to large photometric noise or systematic 
effects caused by non-uniform photometric conditions throughout the survey. We 
use robust variable identification methods to extract 6,520 variable candidates 
from this dataset, resulting in an overall variable fraction of ~ 2.9% at the level 
of ~ 0.05 mag variability. Despite the sparse and uneven time-sampling of the 
light-curve data, we discover 143 periodic variables in total. Due to period am- 
biguity caused by relatively poor phase coverage, we identify a smaller final set 
of 101 periodic variables with well-determined periods and light-curves. Among 
these are 55 RR Lyrae, 30 eclipsing binary candidates, and 16 high amplitude 
Delta Scuti variables. In addition to these objects, we also id entify a sample 



of 2,7 04 variable quasars matched to the SDSS Quasar Catalog (ISchneider et al. 



20071 ). which make up a large fraction of our variable candidates. An additional 
2,403 quasar candidates are tentatively identified and selected by their non-stellar 
colors and variability. A sample of 11,328 point sources that appear to be non- 
variable given the limits of our variability sensitivity is also briefiy discussed. 
Finally, we describe several interesting objects discovered among our eclipsing 
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binary candidates, and illustrate the use of our publicly available light-curve cat- 
alo£ by tracing Galaxy halo sutetructure with our small sample of RR Lyrae 
variables. 

Subject headings: binaries: eclipsing — catalogs — stars: variables: other — 
surveys 



Introduction 



The last fifteen years have yielded a wealth of new knowledge of many types of astro 
nomical obj ects due to the pr evalence of large scale surveys such as the Sloan Digital Sky Sur- 



vey ( SPSS: lYork et al.ll2nnnh and the Two-Micron All Sky Survey (2MASS; ISkrutskie et al. 



20061 ). These surveys cover large areas of the sky, are deep, are sensitive to many different 
types of objects, and most importantly, have uniform data reduction and characterization. 
Powerful data access tools made available to the community allow efficient dissemination 
and data-mining, and allow large populations of objects to be studied all at once. 

The next generation of astronomical surveys will include a powerful new tool: the 
exploration of the time-domain. Many classes of astronomical objects vary over differing 
timescales, so a uniform approach that covers all of these timescales while covering large 
portions of the sky will produce a diverse catalog of variable objects that may be studied 
in much the same manner as the static sky has been studied using SDSS and 2MASS. Two 
such projects in adv anced stages of planning and development are the Pan-STARRS project 
( iKaiser et al.l l2002l: currently undergoing comrn issioning) and the Large Synoptic Survey 
Telescope (LSST: iTysonI l2002l : llvezic et al.ll2008l : expected to start observations in 2015). 



On a smaller scale, the variable sky has been explored by ra any successful projects, in- 
cluding OGLE (toward s the Galactic Bul ge: ludals ki et al.ll2002L MACH O (Galactic Bulg e 



and Magellanic Clouds; lAlcock et al.ll200ll ). and ASAS (all-sky to V = 15. 0:lPoimanski 



Photo metric surveys for t ransiting planet s in recent years, suc h as HAT (IBakos et al 



TrES flAlonso et al.l 120041 ) . SuperWASP flPoUacco et al.ll2006f ). and XO flMcCuUough et al. 



2002). 



2002) 



20061 ) have also resulted in studies of several classes of variable stars. In addition to these 
surveys, specia l ized s earches for supernovae, including the C FHT Supernova Leg acy Survey 
( ISuUivan et al.l l2005l ) and the SDSS-II Supernova Survey (IFrieman et al.l 120081 ) . can pro- 
vide large and uniform datasets suitable for the identification and characterization of many 
different types of variable sources. 



^Available at |http : //shrike . pha ■ j hu . edu/stripe82- variables 
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Although the SDSS is largely a single-epoch survey, it covers some regions of the sky 
multiple times. The most prominent of these is a region on the celestial equator known as 
Stripe 82. This region has been surveyed repeatedly by the SDSS over many years, most 
recently by the SDSS-II Su pernova Survey, and has been studied for variability by several 
authors. ISesar et al.l (120071 ) presented a catalog of 13,051 variable star candidates discovered 
using ob servations of the Stri pe carried out before the advent of the SDSS-II Supernova 
Survey. iBramich et al.l (120081 ) published a light-curve ca talog of Stri p e 82, incorporating 
data from the first two years of the Supernova Survey. iBlake et al.l (120081 ) discovered a 
low-mass eclipsing binary in this region, using l ight- c urves also generated from the first two 
years of the Supernova Survey. iBecker et al.l (120081 ) reported another low mass eclipsing 
binary discove red in the footprint of Str ipe 82, but usi n g cali bration data from 2MASS. 
More recently, IWatkins et al.l (120091 ) and iKowalski et al.l (120091 ) have presented catalogs of 
RR Lyrae and M-dwarf flare stars present in Stripe 82 respectively. 

Here, we construct a light-curve catalog for a magnitude limited sample of point sources 
in SDSS Stripe 82, using data ranging from the first observations of this area of sky in 1998 
to the high cadence observations carried out by the SDSS-II Supernova Survey (2005-2007). 
We use inhomogeneous ensemble differential photometry to remove systematic artifacts from 
these light-curves caused by variable photometric conditions. Point sources are then classified 
by color and by cross-matching to other catalogs of interesting objects in Stripe 82. The 
large number of observation epochs available in our dataset then allows robust identification 
of variable sources, especially periodic variables. 

We first construct a light-curve catalog of variable point sources detected in Stripe 82. 
We then concentrate on three classes of periodic variables: eclipsing binaries, RR Lyrae, and 
Delta Scuti, and present periods and phase-folded light-curves for all such objects identi- 
fied. In addition, we use the color an d variability propertie s of quasi-stellar objects (QSOs) 
matched to the SDSS Quasar Catalog (ISchneider et al.ll2007l ) to identify a new sa mple of vari 
able c andidate QSOs in our light-curve catalog. Finally, we extend the work of llvezic et al. 
( 120071 ) by using a larger number of observation epochs and our robust variable extraction 
methods to identify a sample of objects that appear nonvariahle at the limits of our sensi- 
tivity. 

In this paper (the first of two), we first describe how we extract objects from the de- 
tection catalogs from the SDSS pipeline, and subsequently organize and generate an initial 
light-curve catalog for the RA to 4 h half of Stripe 82 (Section [2]). We describe the in- 
homogenous ensemble differential photometry algorithms in detail (Section 12.31) . and then 
discuss how point source classification and variable extraction are implemented (Sections 
12.41 and l3.ip . Two independent period search algorithms are described and the difficulties 
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we face in searching for periodic variability in our sparse and unevenly sampled dataset are 
also outlined (Section I3.2p . We then apply these variable extraction and period finding al- 
gorithms to objects in the RA to 4 h half of Stripe 82, and describe the general properties 
of variables discovered after processing this initial light-curve catalog (Section 14. 11) . 

We characterize the completeness and efficiency of our variable extraction pipeline by 
carrying out end-to-end simulations of the entire process (Section 14.21) . Our period finding 
algorithms are analyzed in a similar manner and provide estimates for the efficiency and 
completeness of our periodic variable sample. Instructions on how to access our publicly 
available data are then provided (Section l4.3p . We then present and discuss the properties 
of eclipsing binaries, RR Lyrae, and Delta Scuti variables identified in this initial light-curve 
catalog (Section [4.41) . Finally, we discuss samples of candidate QSOs and nonvariable objects 
also identified by our pipeline (Sections 14.51 and 14.61) . The second paper in this series will 
complete our light-curve catalog by adding objects and variables identified in the RA 20 
to h half of Stripe 82, and extend the discussion of their properties in the context of a 
completely processed sample. 

2. Data from SDSS Stripe 82 
2.1. Overview 



The Sloan Digital Sky Survey (SDSS) uses a dedicated 2.5-m telescope located at the 
Apache Point Observatory in New Mexico. During imaging, the telescope points at a fixed 
declination on the meridian and scans the sky using the time-delay and integrate mode that 
clocks the charge across the charge-coupled device arrays (CCDs) at the sidereal rate. The 
record of one particular scan of a ~ 1.3°- wide strip is called a run. Two overlapping strips 
make up one stripe; as a result, each stripe is ~ 2.5° wide. Each run exposes five CCDs in 
five filters, u, g, r, i, and z for approximately 54 seconds per pixel. One run is further broken 
up into fields, which are 13' x 10' each. Detections are fed through the SDSS photometric 
pipeline, are classified by morphological type, and are assigned magnitudes and associated 
uncertainties. Point source objects are well described by fitting a point spread function 
to the detection, resulting in so-called PSF magnitudes, while extended source objects are 
characterized by model m agnitudes. Details of the SDSS photorn e try and classifica t ion a l- 
gorithms rnay be found in iGunn et al.l (120061 ) , iLupton et al.l (120021 ) , IStoughton et al.l (120021 ) , 



Hogg et al.l (l200ll ). and references therein. Observations take place on nights that have see- 



ing better t han 1.7" (FWHM ) , are moonless, a nd show little transparency variation due to 
cloud cover jYork et allboool : Iflogg et allboOlh . 
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Stripe 82 is on the celestial equator, ranging from 20 li to 4 h in right ascension and 
-1.27° to +1.27° in declination, for a total area of about 300 sq. deg. The region is at 
high Galactic latitude, ranging from —30° to —70° in Galactic latitude and between 40° and 
200° in Galactic longitude. The Stripe has been observed many times over the 10 years of 
operation of SDSS and SDSS-II, mostly during the fall (September - December). Figured] 
(left panel) shows the temporal coverage of the Stripe as a function of time over all ten years 
of observation. 

Starting with commissioning runs in 2004, and full operations in 2005, Stripe 82 was 
observed at higher cadence by the SDSS-II Supernova Survey (hereafter, the SN Survey; 



Frieman et al.l 1200811. The nominal time between consecutive observations of a field was 



roughly 2 days (see Figured], right panel). The SN Survey was designed to obtain a uniform 
sample of medium-redshift {z ~ 0.05 — 0.35) Type la supernovae for the purposes of precision 
cosmology. The requirement for improved temporal coverage necessitated a compromise on 
the photometric quality. As a result, the SN Survey observed on nights that would not 
have been designated as photometric for the legacy SDSS. This relaxed photometric quality 
requirement poses a challenge for identification of 'true' variables from the dataset, as we 
will note in Section 12.31 

In our dataset, we have 242 runs from Stripe 82: 69 from observations of the Stripe 
before the SN Survey, 65 for the first season of the SN Survey in 2005, 86 runs during 
the second season in 2006, and 22 runs during the final season in 2007[^. We obtained the 
pre-SN Survey data for Stripe 82 using the Catalog Archive Server (CAS) SQL interfac^ 
to the SDSS database, downloading FITS files with lists of object detections per run along 
with their photometric properties. The data for the SN Survey was reduced by the SDSS 
photometric pipeline and made available on the SDSS Data Archive Servei0 as calibrated 
object catalogs (tsObj FITS binary table files, one per field) and calibrated images (fpC FITS 
image files, one per field). We downloaded the calibrated object catalogs in form of tsObj 
FITS files for runs associated with the SN Survey, put them together with the object catalogs 
for the pre-SN Survey runs obtained earlier, and then processed all 242 runs through a multi- 
stage pipeline designed to extract suitable point sources and search for possible variables. 
An overview of the pipeline is given in Figure [S] while details are discussed below. 



full list of runs used is available at our catalog's website. 
" http : //cas j obs . sdss . org/Cas Jobs/| 




http : //das . sdss . org/www/htinl/imaging/dr-byRuii-74.htinl 
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2.2. Object Extraction and Pre-Ensemble Processing 

The data fields extracted from the SN Survey data are hsted in Table [H We imposed 
a uniform set of quality conditions using status and quality fiags available in the tsObj and 
CAS output FITS fileqj. We first required that a detected object: (1) was not marked as 
saturated (BRIGHT and SATURATED fiags set to 0), (2) was not a deblended child of a 
nearby object (PARENT set to -1, or PARENTID set to 0), (3) had no deblended children 
itself (NCHILD set to 0), (4) had the EDGE fiag set to 0. Second, we required a detected 
object's STATUS fiags included SET, GOOD, OK_RUN, OK_SCANLINE, OK_STRIPE all 
set to 1, in addition to either PRIMARY or SECONDARY set to 1. We explicitly discarded 
all objects that had a status fiag of DUPLICATE set to 1, indicating a duplicate detection 
within the same observation run (in overlapping adjacent fields for example). Third, we 
ensured that an object had valid magnitudes recorded in all five bands by requiring that these 
measured values all be greater than -9999.0. Finally, we imposed a faint PSF magnitude 
limit of z = 21.0. Although we extracted both PSF and fiber magnitudes for each object 
from the catalogs, we only used the PSF magnitudes for all of our subsequent processing 
and variability analysis. 

We then extracted only those objects classified by the SDSS photometric pipeline as 
point sources, by requiring the object TYPE fiag be set to 6 (STAR). This left us with 
57,406,616 point source detections over all runs. We then constructed a canonical match 
template for the entire Stripe using object detections in four high photometric quality SN 
Survey runs: North runs 5610 (MJD 53628) and 6430 (MJD 54011), and South runs 5776 
(MJD 53669) and 6425 (MJD 54010). These runs were chosen because of acceptable seeing 
on these nights (median of 1.7") and full coverage of the Stripe from 20 h to 4 h in right 
ascension. Detections were matched between these runs using a 5.0" radius, and multiple 
detections of the same object discarded. The final match template then contained all point 
sources detected at least once in any of these runs, for a total of 2,000,242 objects. 

We then matched all detections over all runs to the match template using a 5.0" match 
radius. All detections grouped inside this match radius were considered detections of a single 
match template object, and constituted a match bundle. The typical separation between 
neighboring point sources in our match template is ~ 35". This is comfortably larger than 
our chosen match radius, so there is little chance of a detection being mistakenly associated 
with a nearby match template object. We also note that the match radius sets a rough upper 
limit of ~ 0.5" yr~^ for the proper motion of any one match template object, assuming it is 
detected over all 10 years of temporal coverage. 



^See |http : //www . sdss . org/ dr7/products/catalogs/f lags . html 



for details. 
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Figure [3] presents the distribution of the number of detections per match bundle. There 
are 473,759 objects (23.7% of the total) that have fewer than 10 detections in all 242 epochs 
over all 10 years of available photometric data. Of these, ~28.6% have only a single detection, 
~35.4% have extreme color^ most likely due to unreliable photometry, and ~55.8% are 
faintlll. In contrast, only ~0.2% of objects with at least 10 detections have extreme colors, 
while only ~2.3% of the objects with at least 10 detections are classified as faint. We, 
therefore, err on the side of caution and remove all objects with fewer than ten detections from 
any further consideration to keep contamination by unreliable photometry to a minimum. As 
a result, we have 1,526,483 objects with 55,610,252 total detections in our final point source 
light-curve catalog for Stripe 82. We carry this catalog forward to the ensemble photometry 
stage of our processing pipeline, described below. 



2.3. Ensemble Differential Photometry 

As discussed earher, the SN survey sacrifices some photometric quahty to gain temporal 
coverage. Unfortunately, this means that observations are taken under all kinds of photomet- 
ric conditions, including sub-optimal seeing, variable transparency, and at all lunar phases 
except for full moon. We note that the Stripe 82 runs before the SN Survey have more 
stringent constraints on these variables, and are generally of higher photometric quality. 
Even these observations, however, suffer fro m variable photometric conditions on time-scales 



of minutes to hours (see llvezic et al.l 120071 for details). We, therefore, face the challenge of 
reconciling the low cadence but relatively high quality measurements taken before the SN 
survey with measurements of lower quality at much higher cadence taken during the Survey. 

In particular, identification of variable objects becomes problematic when large system- 
atic effects caused by night-to-night sky and seeing variations are present. This is illustrated 
by an example relation between median SDSS light-curve magnitude and light-curve stan- 
dard deviation (Fig HI top-left and top-right panels). This relation is often used to identify 
obvious variable objects; these lie above the general trend. It is difficult to identify true 
variables, however, when the relation itself has a scatter caused by magnitude measurements 
affected by differing photometric conditions. In addition, as seen in the left panel of Figure [5l 
these systematic effects usually appear as dips in light-curves on specific observation dates. 
These dips can be very large and affect measurements in all bands for all objects in a field. 



^Defined as \u - g\ > 3.0 or I5 - r| > 3.0 or \r - i\ > 3.0 or \i ~ z\ > 3.0. 

^Defined as g and r > 22.2, i > 21.3. These are the 95% point source detection repeatabiHty Hmits for 
these bands. See http : //www . sdss . org/dr7| for details. 
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thus artificially increasing the light-curve standard deviation of even completely nonvariable 
objects. 

Such artifacts are usually removed in variability studies by the use of differential pho- 
tometry. A target star is compared to an ensemble of comparison stars, and the measured 
magnitudes normalized to a weighted ensemble magnitude. The resulting differential magni- 
tude for the target object removes any systematic variation in photometric conditions that 
affects all stars being observed on the same night. This technique is powerful and has been 
widely deployed, owing to its computational simplicity, as well as the very high photometric 
precisions that are possible. We, however, cannot employ this simple differential photometry 
method, because we are not guaranteed to have the same comparison stars in the ensemble for 
any one target star over all observations. Furthermore, we must ensure that all comparison 
stars are not intrinsic variables themselves; this is not possible a priori when studying light- 
curves of ~ 1.5 million objects. We turn, theref ore, to a different ial photometry technique 



known as inhomogeneous ensemble photometry (IHoneycuttl 119921 ) . The general method is 



described immediately below; a discussion of its application to our dataset follows thereafter. 

We first assume that most of the stars in a given field containing a target star are 
nonvariable, and that the measured magnitude of any one star z at a time index value of j 
is given by 

m{i,j)=mo{i) + e{j), (1) 

where m{i,j) is the measured magnitude of a star, mo(i) is the 'true' mean magnitude that 
would have been measured in the absence of photometric condition variations, and e{j) 
represents the change in photometric zero-point caused at time index j by these variations 
and is applied to all stars in the field at that time index. We must minimize the quantity (3, 
given by 

Ngtars Nobs 

/5= ^ "^[miij) -moii) - e{j)fw{i,j), (2) 
i=i j=i 

where w{i,j) is the weight associated with a measured magnitude m{i,j) and is given by 

w{ij) = Wi{j)w2ii)w3{i, j)w4{ij). (3) 

The weights wi, W2, and ws are either zero or one, while W4 is the statistical weight based on 
the uncertainty of measurement a[m{i,j)] for m{i,j) and is calculated as l/a'^[m{i, j)]. The 
weight Wi = if the observation at time index j is to be excluded, 'u;2 = if the star i is to 
be excluded, and W3 = if a specific observation of star i at time index j is to be excluded 
from the photometric solution. Once this solution is computed, and we have obtained the 
'true' mean magnitudes mo{i) of all stars in the ensemble, along with the error terms e(j). 
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we can solve for the corrected magnitude M{i,j) for any star i at time index j: 

M{t,j)=m{t,j)-e{j). (4) 
Finally, we can calculate the variance of the 'true' mean magnitudes (7^[mo(i)] using: 

Nobs Ef=i[M{i,j)-mo{z)f wiz, J 



The relation between mo(i) and a[mo{i)] is equivalent to the usual magnitude-cr relation. 
Objects with large values of a[mo(^)] relative to the general trend at that magnitude may be 
considered as possible variables. 

In practice, our inhomogeneous ensemble photometry implementatiorj^ normalizes the 
'true' mean magnitude mf){i) for a star i to the 'true' mean magnitude of the brightest star 
in the ensemble, resulting in effective differential magnitudes M{i,j) for all of the stars. We 
use the weights Wi, W2, and W3 to take into account the presence or absence of comparison 
stars in the ensemble from night to night. To deal with the problem of comparison stars 
being possible variables, we run the photometric solution for a target three times, each time 
removing all comparison stars from the ensemble that appear to be variable. The ensemble 
itself is chosen such that it has at least 10 comparison stars observed within 10 seconds of the 
time of observation of the target star (so within 150 arcseconds of the position of the target 
star, assuming tracking at sidereal rate). The median number of comparison stars is ~ 80 
per target star. This set of conditions ensures that we obtain the best possible differential 
magnitude for the target on each night that it is observed. At this stage, we also remove 
from any further consideration all target stars that do not have comparison star ensembles 
satisfying these conditions, since reliable differential photometry is not possible for these 
objects. In the processing of our RA to 4 h light-curve catalog, this removes only 2, 653 
objects (~ 0.72% of the total) from consideration. 

We show the impact of inhomogeneous ensemble photometry in Figures H] and [51 The 
relation between the median light-curve magnitude and the differential light-curve standard 
deviation (Fig SI bottom-left and bottom-right panel) is seen to be much improved, with 
much of the extrinsic scatter removed. Furthermore, the light-curves of the target object 
and its four closest neighbors (Fig[5l right panel) no longer suffer from the systematic effects 
seen in the left panel. 

The fourth order polynomial fits to the empirical relation between the mean SDSS light- 
curve magnitude (in each band ugriz) and the differential magnitude light-curve standard 



See M. W. Richmond's website: jhttp : //spiff . rit . edu/ensemble 
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deviation a calculated using differential magnitude light-curve data for 365,086 objects with 
at least 10 observations from our initial light-curve catalog (RA to 4 h) are given below: 

a{u) = 105.94 - 24.197m + 2.0621^^ _ 0.07772m3 + 0.00109m^ (6) 

a{g) = 51.441 - 11.337^ + 0.93529^^ - 0.03423c/^ + 0.00047^^ (7) 

a(r) = 48.684 - 10.831r + 0.90210r2 - 0.03334r3 + 0.00046r^ (8) 

a{i) = 54.844 - 12.291i + 1.0315^^ - 0.03842i^ + 0.00054?^ (9) 

a{z) = -69.446 + 15.591^ - 1.0329^^ + 0.04799^^ - 0.00066;z^ (10) 

We generate differential magnitude light-curves in all five bands for all objects in our 
light-curve catalog and search them for variability. The large number of stars in our dataset, 
coupled with the numerous lookups required to build each target star's comparison ensemble, 
make this stage in our pipeline the most computationally intensive, and thus, the slowest. 



2.4. Classification of Point Sources 



We classify the point sources in our final light-curve catalog in two ways. First, we 
cross-match objects between our catalog and existing catalogs of interesting and possibly 
variable objects in the SDSS, especially those present in the Stripe 82 footprint. Second, we 
use several selection algorithms to sort objects into rough classification bins by their SDSS 
colors. In this way, we can select interesting populations of objects to study for variability 
without having to resort to spectra. 

Several authors have already compiled catalogs of interesting classes of point source 
objects discovered during SDSS operations. ISchneider et al.l (120071 ) reported the discovery 
of 77,429 quasars ranging in redshift from 0.08 to 5.41, covering 5,740 sq. deg of sky. Most 
quasars are suspected to be variable in some degree; our catalog of multi-color light-curves 
covering a base-line of nearly ten years may be used to test that assumption. Other catalogs 
that include possible variable objects in clude a catalog o f 9,316 spectroscopically confirmed 
white-dwarfs and 948 hot sub-dwarfs by lEisenstein et al. and the already mentioned 

SDSS Stripe 82 catalog of variable stars by ISesar et al.l (120071 ). containing 13,051 sources 
(hereafter referred to as the SDSS-I Variable Object Catalog). We also cross- match our 
cat alog with th e Strip e 82 Standard Star Catalog (containing ~ 1 million objects) generated 
by llvezic et al.l (120071 ). to confirm the nonvariability of these objects and their suitability 
for use as faint standard stars. We use a 5.0" radius to match between all of the preceding 
catalogs and our catalog of match template objects. 



In addition to classifying objects in this way, we also make use of the excellent color 
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selection made possible by the five photometric bands used in the SDSS. Specifically, we use 
the median SDSS light-curve colors u — g, g — r, r — i, and i — z for each object. We use color 
selection algorithms from three sources: (1) spectroscopic target selection color cuts from the 



SDSS -II Sloan Extension for Galactic Understanding and Expl oration (SEGUE; lYanny et al. 



2009 ) to select different types of stellar objects, (2) colo rs from Sesar et al. ( 20071 ) to identify 
potential RR Lyrae, and (3) color cuts for M-dwarfs in West et al. ( 20081 ). Table [2] presents 
the color classification schemes we use for our catalog. The SDSS photometric pipeline 
reports magnitudes measured for each detection, and the extinction in magnitudes at the 
position of the detection. We use these measurements and calculate the dereddened magni- 
tudes for all objects. These are then used for all color selection cuts, except where noted in 
Table [21 We do not use dereddened colors for brown dwarfs and main-sequence/white-dwarf 
pair selection because these faint objects are likely to be nearby. We note that any one point 
source may be assigned multiple categories based on color; these are all noted in classification 
tags associated with that point source in our catalog. Objects that are not matched to any 
of the catalogs listed above and cannot be classified using color are also noted and assigned 
a classification tag of unknown. 

Results from the catalog cross-matching and color classification schemes above for the 
first release of our catalog are reported in Section HI We now turn to the identification of 
variable objects in our catalog. 



3. Variable Objects in Stripe 82 

3.1. Extraction of Variable Point Sources 

We identify possible variable sources in two steps. The first involves using the results 
from the ensemble photometry stage of our pipeline, specifically the relation between the 
differential magnitude light-curve standard deviation (T[mo(i)] (Equation E]) and the 'true' 
mean differential light-curve magnitude mo{i)- As the ensemble differential photometry 
stage of our pipeline iterates over each target object, it generates a plot of this relation for 
all objects in its associated ensemble, and fits a second-degree polynomial to the empirical 
trend. Outliers more than 2-a away from the general trend are iteratively discarded to make 
the fit more robust. 

This procedure is then repeated to obtain such plots and trends for the u, g, r, i, and 
z bands separately. Examples of these plots are shown in Figure El These are then used to 
tag the target object as a tentative variable if: it lies at least 2-a above the general trend 
of the magnitude- standard deviation relation in (1) u and g bands simultaneously, or (2) 
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g and r bands simultaneously, or (3) r and i bands simultaneously, or (4) i and z bands 
simultaneously. We also tag objects as tentative variables if they lie at least 1-cr above the 
general trend in all three of the r, and z bands simultaneously; this has the potential 
of picking up faint red variables in populations of K and M-dwarfs. This system, however, 
is not robust against false variability caused by large photometric noise for faint objects. 
These objects are likely to lie above the variability threshold due to their small numbers and 
large uncertainties, even in their differential magnitudes, and thus end up being erroneously 
tagged as tentative variables. 

The second step involved in identifying variabl e sources corr ects for this type of false 



variability. We employ the Stetson variability index (IStetsonlll996l ). This is a measure of the 
correlation of simultaneous variability across a pair of bands. The Stetson variability index 
Jai3 for two bands a and [3 observed at the same time is 



k 



JaP — , (11) 

where Wk is the weight assigned to fcth pair of observations. The product of the normalized 
magnitude residuals for the fcth pair of observations is given by 



pa/3 ^ / ^« / ^/3 / Qfc - a \ f Pk- P j ^^2) 

where au and Pk are the magnitudes measured at time index fc, aa,k and a^^^k are the errors 
assoc iated with thes e measurements, a and /3 are the weighted mean magnitudes, iterated 
as in IStetsonI (119961 ) to avoid outliers, and n is the number of observations. The Stetson 



variability index is large when the normalized magnitude residuals are correlated, as 
in the case of real variable sources. A nonvariable source, even one with a large light-curve 
standard deviation, will have uncorrelated magnitude measurements across pairs of bands. 
This will drive the Stetson index for such objects toward zero. 

We calculate the Stetson variability indices J„g, Jgr, Jri, and Jj^ using the differential 
ugriz magnitude light-curves obtained from the ensemble photometry procedure described 
in Section 12. 3[ Taking into account the variability index values for obvious variables in 
the dataset identified by inspection of a large number of differential magnitude light-curves 
(and phased differential magnitude light-curves in the case of periodic variables), we set 
a threshold index value of 0.3 to separate variable and nonvariable objects. Sources are 
therefore tagged as probable variables if: (1) any one of J^g, Jgr, Jri, and Jj^ is greater than 
this threshold value, and (2) Jgr and are both greater than 0.0. Figure [7] shows the 
distributions of the calculated Stetson indices with their standard deviations (dot-dashed 
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lines) and threshold Stetson index value (dashed line) marked. Variable objects are selected 
relatively efficiently using our threshold Stetson index value. 

Figure [S] shows the distribution of Stetson Jg^. with median light-curve r magnitude for 
all 365,086 objects with at least 10 observations in the RA to 4 h light-curve catalog. 
The underlying distribution shows no discernable trend with magnitude, indicating that the 
ensemble photometry routines combined with the Stetson index as an indicator of variability 
are quite robust in variable selection. Furthermore, as we go to fainter magnitudes (beyond 
r ~ 22.0), the distribution of the Stetson index no longer shows an outlying population of 
variables, indicating that these fainter objects have no correlated variability detected and 
thus appropriately have smaller Stetson index values. This points out the effective faint 
magnitude limit of our search for variability (discussed further in Section [4.11 below). We 
quantify the efficiency and robustness of our variable detection routines in Section 142) using 
synthetic light-curve catalogs and associated simulations of variability analysis. 

Probable variables selected by this two part procedure can then be searched for periodic 
variability. If such variability exists, we generate estimates of the period, as described in 
Section 13.21 below. 



3.2. Finding Periods 



We search for periodicity among all objects marked as pr obable variables u sing two inde- 
pendent period-finding methods: the string length method of lDworetskyI (119831). and a varia- 
tion o n the cla ssic Lafle r -Kinm an phase- dispersion minimization algorithm (ILafler fc Kinman 
19651 ) used by [Stetson] (119961 ). Both of these methods attempt to minimize the sum of the 
dispersions of measurements ordered by phase for a test period (the 'string length'), in order 
to produce the 'smoothest' possible light-curve. Perhaps their most important feature is 
that they do not involve the calculation o f sinusoidal components for a light-curve, as in the 
Lomb-Scargle periodogram (IScarglelll982l ). and thus do not pre-suppose such a shape when 
searching for periods. This makes them invaluable for studying a broad range of periodic 
variable types, as expected in our dataset. These two algorithms also do not involve binning 
consecutive li ght- curve measurements, as se en in more sophisticated period finding methods 
such as AoV (ISchwarzenberg-Czernyl Il989l ). and BLS (IKovacs et al.ll2002l ). BLS and AoV 
do not work well on datasets such as ours, which have light-curves with a very small number 
of unevenly distributed and sparse time-series measurements over a long baseline. 



The Dworetsky string length Ldw calculated for a test period and a series of phase- 
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ordered photometric measurements mi to is given by the relation 

N-l 



Ldu, = J2 V'K^^^^Ki^TiiiT^^^ + ^K-m'^)2 + (01 -0^ + 1)2, (13) 
1=1 

where (pi is the phase of observation i, and m' are modified magnitudes used to assign similar 
weights to the time and magnitude measurements, and are given by the relation 

m- = {rrii - rrimin) / {2.0{mraax - rrimin)) - 0.25, (14) 

where rrimin and rrimax are the minimum and maximum magnitude measurements in the 
timeseries, respectively. 

The Stetson algorithm incorporates the weighting of individual light-curve measure- 
ments by their respective errors, thus, providing a more robust means of calculation of the 
string length. The Stetson string length Lst calculated for a test period and a series of 
phase-ordered photometric measurements mi to mjy is given by the relation 



Ei=l^(«.« + 1) 



(15) 



where w are the weights assigned to magnitudes m, which in turn have measurement errors 
a. The expression w{i,i + 1) is given by the relation 

where cTj is the measurement error associated with magnitude measurement m^. 

We calculate L^w and Lgt using griz differential magnitude light-curves for all objects 
tagged as probable variables using a test period interval of 0.1 to 100.0 days (see Figure [9]for 
an example) and a frequency step-size of 0.0001 days~^. The u-band differential magnitude 
light-curves have poor signal-to-noise for many of the objects in our dataset, and are therefore 
not used for string length calculations. Objects with light-curve variations on longer time- 
scales (i.e. years) do exist in our dataset, but we do not attempt to fit periods to these, mainly 
due to extremely poor phase coverage (typically 40-60 unevenly sampled measurements over 
10 years). If, upon inspection of its phased light-curve, an object appears to be variable on 
a time-scale shorter than 0.1 day, we rerun the period finding routines using a test period 
interval of 0.01 to 10.0 days. 

For each object, we retain the test periods that result in the 20 shortest string lengths 
for both Ldw and L^t independently, and then phase-fold the hght-curves using these periods. 
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The most likely period in each case is the one with the shortest string length. We require that 
the two period finding methods agree on the most likely period before accepting an object 
as a likely periodic variable. Further, we make use of the available multi-band photometric 
data by requiring that the most likely period reported by the string length algorithms be the 
same for each of the gri light-curves in the case of non-M dwarf objects, and for each of the 
riz light-curves for redder objects such as M-dwarfs. We then attempt to refine the period by 
rerunning the string length algorithms in a small period interval (typically 0.1 day) centered 
around the already determined most likely period. Finally, we visually inspect each phased 
light-curve obtained using this refined most likely period to assess its credibility, and place 
the object into one of three classification bins based on its light-curve shape: (1) eclipsing or 
ellipsoidal binary, (2) sinusoidal variable with an asymmetric light-curve, (3) and sinusoidal 
variable with a symmetric light-curve. 

As mentioned earlier, our dataset contains light-curves with a small number of measure- 
ments scattered unevenly over large baselines. A periodogram of the typical spectral window 
function for our data is shown in Figure [TOl There is significant power in the peaks around 
1.0, ~ 7.0, and ~ 30.0 days and their aliases, which are related to the three main duty cycles 
present in the SN Survey dataset. Fortunately, string length period-finding algorithms are 
relatively insensitive to these cycles (as evidenced by Figure[9]), but objects with periods close 
to the periods of these cycles suffer from severe aliasing, making it difficult to distinguish 
between the many likely periods. We also face the challenge of insufficient sampling near 
points of maximum variation in light-curves; for example, eclipsing binary candidates that 
do not have light-curve points that sample primary and secondary eclipses will suffer from 
ambiguity in period determination. Finally, objects that have large photometric noise even 
in their differential magnitude light-curves, such as faint red stars, will return many different 
string lengths that do not correspond to any 'smooth' phased light-curve, thus making it 
impossible to determine their periods, even if evidence of periodicity is apparent in these 
objects' unphased light-curves. 



4. The RA h to 4 h Light-curve Catalog 

4.1. General Properties and Classification 

The first release of our catalog covers the right ascension range to 4 h. In this region, 
there are 495,797 objects extracted from the SDSS dataset as described in Section |2] above. 
We restrict our attention to only those objects with at least 10 observations, and process 
these through our ensemble photometry pipeline discussed in Section 12. 3[ The resulting 
light-curve catalog contains 365,086 point sources. We then remove objects that have large 
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errors in color: a{u — g) > 0.2, a{g — r) > 0.2, cr(r — i) > 0.2, ar id a(i — z) > 0.2. This 
pruning leaves 228 ,056 objects. The SDSS imaging pipeline (see iLupton et al.l |2002| and 
Lupton et al.ll200ll ) efficiently separates stars and galaxies up to r = 21.5 mag. To minimize 
the contamination fraction of galaxies in our sample, but at the same time remain sensitive 
to faint red variables, we impose a faint magnitude limit of r = 22.0 mag on the objects in 
our light-curve catalog. This leaves us with a final catalog of 221,842 objects to consider for 
variability analysis. 

We then classify the objects in this light-curve catalog by color-selection and cross- 
matching against other catalogs as discussed in Section 12.41 The results of this process are 
shown in Table [31 M-dwarfs make up the largest fraction of point sources by number; partly 
due to the initial single magnitude cutoff imposed at z = 21.0 mag (see Section [2]), and partly 
because of the intrinsic frequency of these low mass stars in the Galaxy. These objects are, 
therefore, overrepresented in our catalog to a signifi cant degree. We also cross-match 4,196 
objects classified as QSOs by lSchneider et al.l (120071). 419 spectroscopically confirmed white 



dwarfs, and 15 hot subdwarfs from the catalog of lEisenstein et al.l (120061 ). 165,591 point 
sources in ou r catalog are succes sfully cross-matched to objects classified as standard star 



candidates in Ivezic et al. 



variable candidates by ISesar et al.l (120071 ). 



(120071). Final ly, we recover 3,972 objects classified as Stripe 82 



Our variable extraction routines produce 6,860 (~ 3.1% of the total objects) tentative 
variable candidates, and a final list of 6,520 (~ 2.9% of the total objects) probable variable 
candidates in our dataset. Figure [TT] (left panel) presents the variable fraction as a function 
of SDSS r magnitude. This fraction rises to a maximum of ~4.5% at the level of ~ 0.05 
mag for the bin 20.0 < r < 21.0, which contains a significant number of point sources. 
The fall in the variability fraction observed in the bin 21.0 < r < 22.0 reflects the effect 
of increasing photometric noise on the determination of variability. In this magnitude bin, 
correlated variability across several bands is difficult to detect, thus the Stetson index trends 
towards zero, leading to a corresponding decrease in the variability fraction. Overall, we find 
that ~ 2.9% of the point sources in our dataset show variability at the level of ~ 0.05 mag, 
given a median r magnitude of ~ 19.9. The right panel of Figure [TT] shows the trends of 
the differential r light-curve standard deviation with median SDSS r light-curve magnitude 
for nonvariable (solid line) and probable variable (dashed line) objects respectively. The 
large difference between the two for all magnitude bins indicates that variables are robustly 
selected by our methods over this range of magnitudes. 



Figure [T2] shows a g — r/u — g color-color diagram for all point sources (left panel) and 
just the variable candidate point sources (right panel) in our final light-curve catalog. We 
can place all detected sources into six broad classification regions as depicted in this figure. 
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Region A is where white dwarfs lie in g — r /u — g color space, owing to their blue colors. Low- 
redshift quasars cluster in region B, also due to their blue colors. Region C contains A main 
sequence stars, blue horizontal branch stars, as well as blue stragglers. Region D includes 
mostly quasars with high redshifts {z > 3.0) as identified by cross-matching with the SDSS 
Quasar Catalog. Region E is the main stellar locus, with stars getting redder in color space 
as their mass decreases. This results in a progression of spectral types from approximately 
F to early M towards the top-right corner of the color-color diagram. The lowest-mass stars 
(the late M dwarfs), brown dwarfs, and other faint red objects with unreliable u — g colors 
are found in region F. Their u ~ g colors are unreliable due to the large photometric error 
in these bluer bands caused by the very red spectral energy distributions of these objects. 

The right panel of Figure \n\ shows the g — r /u — g color distribution of 6,520 probable 
variable candidates. The two dominant classes of variables appear to be the QSOs in region 
B, and a significant amount of faint red objects present along the stellar locus in region E. 
We note here that the majority of variable objects classified as 'unknown' in Table [3] are 
located within the color space defined by region B, and thus are themselves likely to be 
QSOs. We discuss these further in Section IT5l 

We identify 143 periodic variables among the 6,520 probable variable candidates ex- 
tracted from the catalog. Of these, 12 appear to show periodic variability upon inspection 
of their unphased light-curves, but have extremely poor phase coverage or too much photo- 
metric noise to allow any period to be determined and assigned to the object. Thirty more 
objects, most of which are mid to late M-dwarfs, show eclipse-like periodic variability, but 
have ambiguous periods or phased light-curves that do not show convincing periodicity upon 
closer inspectioio. Figure [13] shows light-curves of a sample of 6 such problematic objects. 
We remove these 30 objects as well as the 12 that have no assigned periods from any further 
consideration. 

We are finally left with 101 periodic variables identified from our light-curve catalog. 
Table H] lists these objects as classified by the shape of their light-curves. The positions of 
these periodic variables inu — g/g — r color space are shown in Figure [HI All of these objects 
are found on the stellar locus, and range in spectral type from A to late M. The RR Lyrae 
and Delta Scuti candidates are clustered in region C, where the blue horizontal branch stars 
may be found, as expected. In contrast, we find eclipsing variables all along the stellar locus. 
No periodic variables are found among the cross-matched white dwarfs or QSOs. 



A list of these is available from 



http : / / shrike . pha. jhu . edu/ stripe82- variables 
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4.2. Variable Detection and Period Recovery Efficiency 

We test our general variable identification process by constructing a synthetic light- 
curve catalog representative of the objects present in our dataset. Running this light-curve 
catalog through our ensemble photometry, variable extraction, and period finding routines 
then provides upper limits for the efficiency and reliability of our detection methodology. We 
first draw median light-curve magnitudes, the associated standard deviations, the number, 
and dates (MJD) of observations for each object from probability distributions matching 
the observed distributions of these quantities. Light-curve measurements are then perturbed 
using Gaussian noise appropriate for the objects' assigned magnitudes obtained from the 
observed magnitude-cr relation (for example. Figure H]). A catalog of 46,000 synthetic objects 
covering an area of 5 degrees in right ascension and 2.5 degrees in declination is constructed 
using these distributions and serves as a 'background' for ensemble photometry of two classes 
of objects: sinusoidal variables and completely nonvariable objects. We insert 2,000 objects 
of each type into the synthetic light-curve catalog, making for a total of 50,000 objects that 
are then processed through our pipeline exactly in the same manner as real objects from the 
Stripe 82 dataset. 

The 2,000 artificial sinusoidal variables are inserted into the light-curve catalog at uni- 
form random spatial coordinates, thus distributing them evenly over the 12.5 square degree 
area under consideration. Periods for these objects are chosen from a uniform distribution 
in logP (where P is the period in days) ranging from 0.1 to 100.0 days. Amplitudes are 
assigned to these objects from a uniform random distribution ranging from 0.01 to 0.50 mag- 
nitudes. We simplify matters by assuming the same variability amplitude for each band. The 
epoch of minimum light for each periodic variable is also chosen from a uniform random dis- 
tribution of observation dates present in the dataset. We then generate sinusoidally variable 
magnitude measurements distributed over the assigned observation dates using the periods, 
amplitudes, and epochs of minimum light thus chosen. Finally, each object's assigned mag- 
nitudes are perturbed by an amount chosen from a Gaussian distribution centered at 0.0 
and with a standard deviation equal to the assigned artificial light-curve standard deviation, 
thus introducing 'noise' to the light-curves of each object. 

Objects are required to have at least 10 observations in their timeseries to be eligible 
for ensemble photometry; 1,587 of the inserted variables survive this cut. Objects with color 
errors greater than 0.2 and magnitudes fainter than r = 22.0 mag are then discarded as in our 
usual processing. We then run ensemble photometry and variable object extraction routines 
on the remaining objects. This leaves us with 1,393 objects to extract variables from. The 
results of these routines are summarized in Table [51 The pipeline identifies 1,101 objects as 
tentative variables] 1,091 of these are eventually tagged as probable variables, resulting in an 
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overall (ideal case) variable recovery rate of ~ 78.3%. The variable extraction procedures are 
most efficient at brighter magnitudes (reaching maximum efficiency near r ~ 18.5). They 
become progressively less so with fainter magnitudes, due to increasing photometric noise 
that makes it more difficult to identify variability. 

The final step in characterizing the efficiency of variable extraction is to test the re- 
covery of the assigned periods for these sinusoidal variables. We carry this out by running 
all detected synthetic probable variables through our period finding routines as described 
in Section 13. 2[ We expect the recovery rate here to be quite low, given the limited phase 
coverage and small number of observations of each object. A periodic variable is considered 
to be recovered with the correct period if the difference between the recovered period and 
the original assigned period is less than 0.001 days. Larger differences introduce signifi- 
cant phasing errors when one attempts to construct phased light-curves for these objects. 
Overall, 600 out of 1,091 artificial sinusoidal probable variables have their periods recovered 
successfully, resulting in a recovery efficiency of ~ 54.9 %. Figure [15] shows how this recovery 
fraction behaves as a function of median r magnitude (top left), number of observations (top 
right), variability amplitude (bottom left), and input period (bottom right). The period of 
the variable object and the number of observations determine the phase coverage, and affect 
the period recovery fraction the most. Long period variables have poor phase coverage and 
a correspondingly small recovery fraction. The fraction drops with increasing magnitude, as 
expected, due to the increasing noise in the light-curves that makes it difficult to phase them 
correctly. The variable recovery fraction generally increases with the number of observations, 
except for the last few bins where the trend becomes ambiguous due to the small number of 
objects in these bins. 

Using results from our simulations, we can also address the question of how our period 
finding algorithms fail when they do. Figure [TBI shows the relation between the assigned input 
periods and the actual recovered periods for all 1,091 artificial periodic variables identified 
by our pipeline. The curves depict the three most common types of aliases for our periods 
present in this dataset. We see that in most cases, the recovered period corresponds closely 
to the input period, but there are significant numbers of cases where the period finding 
algorithms latch on to these aliases in lieu of the actual input period. This is mostly a 
function of the phase coverage on each variable object; the better the balance between an 
object's period and the number of its observations, the more likely it is that we will recover 
the correct period. 

We see this effect most clearly in our sample of actual eclipsing binary candidates identi- 
fied from Stripe 82 data, where we discard nearly half of our initial sample due to ambiguous 
period recovery before settling on a final sample of 30 objects. The rapid changes in these 
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objects' light-curves near eclipse and the very small amount of time each object spends in 
eclipse cause significant difficulties for our period finding algorithms. A photometric followup 
campaign to better characterize these objects would therefore be most successful if it used 
ephemerides generated using the periods reported here as well as several harmonics of these 
periods. 

Finally, we inject 2,000 artificial completely nonvariable objects into the synthetic light- 
curve catalog. These objects only have Gaussian noise added to their light-curves, but 
otherwise show no trends over time. These are used to test the false positive variable identi- 
fication rate of our pipeline. 1,578 such objects survive the initial ensemble pipeline cut on 
number of observations, while 1,399 of these survive our additional cuts on the color errors 
and imposed r magnitude limit. Running the pipeline on these remaining objects yields 12 
objects tagged as tentative variables; only 2 of these are subsequently identified as probable 
variables. Neither of these remaining objects appear to have to have any periodicity, as 
expected. On average, therefore, ~ 0.15% of all artificial nonvariable objects inserted into 
the catalog are expected to be misidentified as actual variables. The actual false positive 
rate will be higher than this number, but we believe our pipeline separates variables from 
nonvariables efficiently, based on the evidence from these simulations. We expect the con- 
tamination fraction of such falsely tagged nonvariable objects in our variable sample to be, 
at worst, near 10%. 



4.3. Obtaining the Light-curve Catalogs 

A complete set of light-curves for the 6,520 variable candidates, lists of the different 
types of periodic variables identiffed here, catalogs of the objects discussed in the following 
sections, and ffnally, summary ffles describing all 221,842 point sources in this dataset, are 
available at the following website: 

http:/ /shrike.pha.jhu.edu/stripe82- variables/ 

All object catalogs are in standard FITS binary table format. Light-curves are provided 
as individual FITS tables, CSV, and PDF ffles associated with each variable object, and 
include various related diagnostics. Details on the format and content of these data ffles are 
available at the website mentioned above. Finally, the software code and programs used to 
construct our light-curve catalog and carry out ensemble photometry are also provided. 
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4.4. Periodic Variables 

4.4- i- Eclipsing and Ellipsoidal Binary Candidates 

We find 30 eclipsing binary candidates in our dataset. The period distribution for all 
of these objects is shown in the form of a histogram in Figure [T71 Example light-curves 
for these objects are presented in Figure [181 and a listing of all eclipsing binary candidates 
found by our pipeline is presented in Table El The comments column in that table indicates 
the type of object as classified by its color, as well as noteworthy features of its light-curve. 
Objects with periods greater than 1.0 day are relatively scarce. This is primarily due to 
the sampling cadence of the survey. Objects with longer periods spend little time in either 
primary or secondary eclipse, and sparse, uneven sampling of their light-curves, coupled 
with a small number of observations makes it difficult to obtain unambiguous periods. In 
contrast, objects with small periods are more suited to our sampling cadence; these are more 
likely to be in eclipse at any one time, and thus provide a better estimate of the period. The 
abundance of short period objects, therefore, is a result of observational bias, and cannot be 
ascribed to physical reasons. 

The binary candidates appear to be rather diverse in their nature. We find at least 
five W UMa type contact binary systems; two good examples of these are objects MB5010 
(SDSS J035138.50-003924.5) and MB6467 (SDSS J031021. 22+001453.9) in Figure [H Ob- 
ject MB23368 (see FigurelH SDSS J025953.33-004400.3) is a very short period late M-dwarf 
binary candidate, which shows an interesting out of eclipse variation in its light-curves, most 
prominently in the g and r band. This may be evidence of a star spot rotating in and out 
of view, but a more densely sampled light-curve for this object is required for confirmation. 

Other interesting binary candidates include objects MB42018 (SDSS J032515.05-010239.7) 
and MB14172 (SDSS J024255.78-001551.5). MB42018 is a mid M-dwarf object that appears 
to have a short period and rather shallow eclipses (~ 0.3 mag for the primary, ~ 0.1 mag for 
the secondary), indicating a possible low mass companion. We note, however, that there is 
a third star present near the eclipsing binary candidate (within 1.0") and this may instead 
indicate that the shallow eclipses are caused by blended light from this system. At the other 
end of the size scale, MB14172 is a blue star with a long period and a deep primary eclipse 
and a much shallower secondary eclipse. It is tagged as a low-metallicity object and possibly 
a giant star by our pipeline, perhaps indicating a giant-dwarf binary system. Finally, we 
point out the well-sampled light-curve of the object MB8125 (SDSS J030753.52+005013.0, 
a short period K dwarf eclipsing binary candidate) as an example of the accuracy in period 
finding possible even with a relatively small number of observations (115 in this case) over 
the ten year length of the survey. Light-curves of all three objects are shown in Figure [181 



-22- 



We note here that, in addition to the newly discovered binaries described above, we 
have also recovered two confirmed low mass eclips ing; binaries presen t in the footprint of 



Stripe 82: the objects 2MASS J01542930+0053266 ( Becker et a: 



Blake et al 



20081 : MO+Ml), and SDSS 



20081 : M4). We were, how- 



J031823.88-010018.4 (also known as SDSS-MEB-1; 
ever, unable to recover the periods of these objects as reported in the literature, because of 
limited phase coverage and a small number of observations, but our variable identification 
method was sufficient to pick out these objects as obvious variable candidates. These two 
objects are listed in Table El along with our own eclipsing binary candidates. 

The poor phase coverage of most of our light-curve catalog, especially for our eclipsing 
binary candidates, precludes detailed analysis of these objects at this point. Dedicated 
photometric and spectroscopic foUowup will be required to draw meaningful conclusions 
about their physical properties as derived from their light-curves. We will discuss results of 
such followup for some of these binary candidates in a future paper. 



4-4-2- RR Lyrae and Delta Scuti Variables 

We find 71 sinusoidal variables in our dataset. These are distinguished from the eclips- 
ing binaries discussed previously by inspection of their phased differential magnitude light- 
curves. These variables can be further classified into three broad types based on light-curve 
shape, period, and amplitude; the RRab RR Lyrae, the RRc RR Lyrae, and high amplitude 
Delta Scuti variables. We fit Fourier components to the light-curves of all 71 sinusoidal 
variables to robustly classify them as belonging to one of the three sinusoidal variable types, 
thus 

N 

/(t) = v4o + ^Afcsin[fccj(t-to) + 0fc] (17) 

k=l 

where Aq is the mean magnitude, Ak is the amplitude of the kth Fourier term, u is the 
angular frequency and is given hj u = 271 /P for a period P, is the epoch, and 0^ is the 
phase for the /cth Fourier term. We restrict the maximum order of the fit iV to 3, because 
of the low number of observations per object, which gives us poor phase coverage. We then 
calculate the Fourier parameters R21 and 02i, given by the respective relations 

R21 = 4^ and 021 = 02 - 201. (18) 
The calculation of these two parameter s allows us t o quantitatively distinguish between 



RRab, RRc, and Delta Scuti variables (jPorettil 1200 ll ). We further require that an object 



be fit successfully using Fourier components to be classified as either one of these three 
sinusoidal variable types and be accepted as such into our periodic variable catalog. Figure 
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shows plots of i?2i against period (left panel) and (j)2i against period. The three classes 
of sinusoidal variable separate cleanly in both diagrams. 

There are 55 RR Lyrae variables identified during this process. Of these, 36 are of the 
subtype RRab, and 19 are of the subtype RRc. Table [7] gives a listing of these objects, their 
median SDSS light-curve magnitudes, and their periods. The periods for the RRab variables 
range between ~ 0.476 and ~ 0.709 days, with a median period of ~ 0.608 days. These 
objects show the expected trend of decreasing amplitude of variation and more symmetric 
light-curves with increasing period. This trend is apparent in Figure [201 which presents gri 
differential magnitude light-curves of a sample of RRab variables identified in this dataset. 
The 19 identified RRc variables have periods ranging from ~ 0.201 to ~ 0.409 days, with 
a median period of ~ 0.301 days. Figure [21] presents example gri differential magnitude 
light-curves for these objects. 

Although the selection of RR Lyrae by light-curve shape and period is more robust than 
selection by color alone, the completeness of our sample of such objects suffers due to poor 
phase coverage. The 55 objects selected for our sample are a meagre fraction of the 125 RR 
Lyrae probable variable candidates selected by color, and an even smaller fraction of the 
1,664 RR Lyrae selected by color alone without regard to their variability. We quantify the 
completeness and efficiency of our selection process by carrying out simulations similar to 
those discussed in Section [T2l 

We generate 2,000 RRab and 2,000 RRc synthetic objects using the magnitude, magni- 
tude error, number of observations, and observation date distributions present in our dataset. 
Colors are assi gned to thes e obje cts using the RR Lyrae candidate color cuts described in 



Section [2^ and llvezic et al.l (120071 ). Periods are assigned to the RRab objects from a uniform 
distribution of periods between 0.5 and 0.8 days, and to the RRc objects from a uniform 
distribution of periods between 0.2 and 0.5 days. Similarly, we assign variability amplitudes 
from uniform distributions of amplitudes between 0.4 and 1.0 mag for RRab, and 0.1 to 0.5 
mag for RRc objects respectively. Light-curves for these objects are generated and they are 
then inserted into separate catalogs of 48,000 nonvariable synthetic objects each, and then 
run through the ensemble photometry and variable extraction routines. We test the recovery 
of these artificial variables in two steps: (1) extracting them as tentative and subsequently 
as probable variables by using our variable identification methods, and (2) recovering their 
assigned periods by using our period finding algorithms. 

Of the 2,000 synthetic RRab variables inserted into the catalog, 1,533 survive the ensem- 
ble process and subsequent conditions on color errors and magnitude limit. Of these objects, 
1,517 (~ 98.9%) are successfully identified first as tentative variables and then as probable 
variables. Overall, 1,312 of these variables are recovered with the correct periods, making 
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for an average recovery fraction of ~ 86.4%. Figure [22] shows the trend of the period recov- 
ery fraction with g magnitude, number of observations, amphtude, and period for synthetic 
RRab variables (sohd hues) as well as RRc variables (dashed lines). Period recovery quickly 
becomes more successful as the number of observations for a given object are increased, 
indicating the important role of good phase coverage. The recovery rate appears to be rel- 
atively insensitive to object magnitude, discounting the objects in the brightest magnitude 
bin, which are much fewer in number compared to those in the other magnitude bins. The 
period recovery fraction also appears to be insensitive to the variability amplitude; RRab 
variables have amplitudes that are large compared to our survey's sensitivity to variability. 
The weak trend with period is reflective of the small period range being explored. 

In contrast to the high variable recovery rate of synthetic RRab variables, the synthetic 
RRc variables suffer due to their small variability amplitudes. Of the 2,000 such objects 
inserted into the catalog, 1,543 survive to the variability analysis stage after cuts on color 
error and magnitude limits. Only 1,132 of these objects are picked up as tentative variables, 
and subsequently as probable variables, resulting in a relatively poor variable recovery rate of 
73.4%. The variable recovery rate shows a significant decreasing trend with fainter magnitude 
after peaking in the magnitude bin 19.0 < g < 20.0, similar to the results of general variable 
recovery simulations presented in Section 14. 2[ The period recovery rate for synthetic RRc 
variables is similar to that for synthetic RRab variables: 934 of 1,132 of these (~ 82.5%) 
are recovered with the correct periods. Figure [22] (dashed lines) shows how this rate is 
related to various input parameters. Once again, the most important factor is the number 
of observations, and by extension, the phase coverage for our variable objects. 

Given the results of the variable and period recovery simulations above, we estimate 
a completeness of no better than 65% for our RR Lyrae sample. Although this is low, 
the objects that are discovered are very likely to be real RRab and RRc variables due to 
successful classification by light-curve shape, period, and Fourier decomposition. We can, 
therefore, attempt to trace halo substructure using these objects for illustrative purposes. 
Assuming an absolute magnitude My = 0.7 for RR Lyrae stars, we calculate the distances 
to the 55 such objects in our catalog. Figure [23] shows a plot of the distribution of these RR 
Lyrae (36 RRab, and 19 RRc) as a function of the distance and the right ascension. The 
large clump near 25 kpc and ranging from 30° to 40° in right ascension is associated with part 



of the Sagittarius dwarf stream (S167-54-21.5 in [Newberg et al.[[2002l ). Other clumps can be 
seen in the spatial distribution plot, and these may be associated with halo substructure as 
well. A detailed examination of RR Lyrae in Stripe 82 and how they relate to Milky Way 
substructure at large distances is, however, beyond the scope of this work, largely due to 
the small number o f RR Lyrae in our sample. Excell e nt tre atments of these topics may be 



instead be found in [de Led (120081 ) and [Watkins et al.l (l2009l ) 
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The final class of sinusoidal variables we discuss are the high amplitude Delta Scuti 
(HADS) variables. These have short periods (< 0.3 days) and amplitudes (up to 0.5 mag) 
greater than those of the usual Delta Scuti stars. We find 16 candidates for such objects 
in our dataset. Table [H] gives a listing of these, along with median light-curve magnitudes 
and periods. Figure [2l] shows example light-curves for 6 of these objects. The median 
period for these objects is ~ 0.063 days, with a minimum period of ~ 0.050 days, and 
a maximum period of ~ 0.093 days. Longer period HADS variables can potentially be 
confused with short period RRc RR Lyrae in both period and color space (see Figure [T^ . 
however, Fourier decomposition provides a robust mechanism of distinguishing between the 
two kinds of variables, as seen in Figure [191 



4.5. Variable Quasars 



The vast majority of quasars are expected to be variable on both short and long 



timescales (see iHawkinsI |2002| . Ide Vries et al.l |2003| . and references therein) . The short term 



variability of these objects is likely to be associated with eruptive events and manifests as 
large-amplitude variations in light-curves (for example, BL Lac type objects). In contrast, 
the long term variability of quasars is dominated by small-amplitude variability, often result- 
ing in increasing or decreasing trends in brightness on multi-year timescales. It is possible to 
use this intrinsic variability to distinguish between objects on the stellar locus and quasars 
(ISesar et al.l 120071 . and references therein). 



Our light-curve c atalog contains 4,196 quasars matched to the SDSS Quasar Catalog 
(jSchneider et al.ll2007l ). These objects mostly fall within regions B and D of the g — r/u — g 
color-color diagram (see Figure [T21 left panel). The objects in the quasar catalog were 
selected using SDSS photometry and spectra and have measured redshifts available. We 
identify 2,704 objects among these matched quasars as probable variables tagged by our 
pipeline (see Figure [25] for example light-curves), resulting in a variable fraction of ~ 0.64. 
Furthermore, these variable quasars make up a sizeable fraction, ~ 0.41, of all identified 
candidate variable sources. 

New QSO candidates in our sample may be identified by taking advantage of their 
intrinsic variability and non-stellar colors. As noted in Section 14. H there are 8,463 objects 
tagged as 'unknown' in our catalog due to the lack of any corresponding objects in any other 
catalogs, as well as no classification posssible by SDSS color alone. Of these objects, 1,102 are 
tagged as probable variables. A large fraction of these probable variable 'unknown' objects 
have non-stellar colors and may be found in region A of Figure [12] (left panel). We identify 
new quasar candidates by requiring that they be tagged as probable variables and satisfy 
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either one of two color cuts: low- 2; quasars using —0.3 < u — g < 0.7 and —0.3 < g — r < 0.5, 
and high-z quasars using u~g > 1.4 and u — g > 1.6{g — r) + 1.3A. These color and variability 
cuts result in the identification of 2,403 QSO candidates. 

Figure [26] (left panel) shows the relation between Stetson indices Jug and for variable 
QSOs matched to the SDSS Quasar Catalog a nd stellar locus objects. Q uasars appear to be 



more strongly variable at shorter wavelengths (jVanden Berk et al.ll2004j ). and this tendency 
shows up in slopes of the Jug vs J^ trends plotted in the figure. Variable stellar locus objects 
show most of their variability at longer wavelengths, and thus have larger relative Jgr, and Jri 
values. The addition of 2,403 variable QSO candidates identified by our color and variability 
cuts to the plot (Figure [261 right panel) does not appreciably change the slopes of the two 
distributions, indicating that these objects largely follow the trend for matched QSOs from 
the SDSS Quasar Catalog. 

Despite these candidate objects matching the colors of, as well having similar variability 
properties to the already known QSOs, confirmation of their quasar nature would require an 
extensive spectroscopic campaign. Lists of all objects matched to the SDSS Quasar Catalog, 
as well as those objects identified as candidate QSOs here are both available from the website 
mentioned in Section WTLl above. 



4.6. Nonvariable Objects 



There are 214,982 objects in the final light-curve catalog (the r < 22.0 sample) that 
fail to meet the selection criteria for tentative variables. Not all of these will be actual 
nonvariables. We can, however, use the tools developed here to define a sample of point 
sources that show no variability at the limits of our sensitivity. Probable nonvariables are 
thus selected if they: (1) are not tagged as tentative v ariables by our p ipeline, (2) do not 
match to any objects in the SDSS-I variable star cata log (ISesar et al.|[2007[ ). (3) have a match 
in the SDSS standard star catalog ([ivezic et al.[[2007l) . (4) do not match to any objects in the 
SDSS quasar catalog ( Schneider et al. 2007I ). and finally, (5) have Stetson variability indices 
Jug, Jgr, Jri, and Jiz all less than 0.05. These conditions select 19,704 objects. A further 
refinement can be made by demanding that these objects have at least 50 detections each, 
ensuring that the Stetson index selection for nonvariability is valid over at least three years 
of observations. With this final criterion in place, the probable nonvariable sample consists 
of 11,328 point sources. 



Figure [27] shows a. g — r /u — g color-color diagram of these objects. Nearly all of them 
are found on the stellar locus and are well-distributed in color from blue A stars to red M 
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stars. These objects form a subset of the sources found in the SDSS standard star catalog 
and have been selected as nonvariables using a much longer timeseries baseline. Additional 
high cadence photometric monitoring would be required to rule out variability below the 
level of what we can detect with our methods {a < 0.05 mag) before these are actually 
deemed suitable for use in a canonical standard star catalog. A list of these stars is available 
from this paper's accompanying website. 

5. Conclusions 

We have constructed a light-curve catalog of 221,842 point sources in the RA to 4 
h half of Stripe 82, and identified 6,520 candidate variables. Of these, 2,704 turned out to 
be already identified quasars, while another 2,403 were classified as QSO candidates due 
to their colors and variability properties. We found 101 periodic variables in this dataset, 
including 30 candidate eclipsing binary systems, 55 RR Lyrae and 16 high amplitude Delta 
Scuti candidates. We also identified a sample of 11,328 point sources that do not appear to 
be variable, based on observations over a long time baseline and rejection from our variable 
extraction algorithms. 

The use of inhomogeneous ensemble differential photometry and the Stetson variability 
index was crucial in identifying possible variables among the objects in the dataset, while 
removing many sources that appeared to be falsely variable. The poor phase coverage of our 
light-curves presented difficulties in extracting periodic variables from the candidate variable 
sources. Binless phase dispersion minimization methods, such as the Dworetsky and Stetson 
string length algorithms worked well on our sparse and unevenly sampled data, and had the 
added advantage of being unbiased with respect to the kind of variability being probed. 

We have made public our Stripe 82 variable object light-curve catalog, along with the 
software implementation of our ensemble photometry pipeline. In an upcoming paper, we 
will address variability in the remaining half of Stripe 82 (RA 20 to h), and present our 
final catalogs for periodic and other types of variables present in this dataset. Unfortunately, 
the relative faintness of many interesting periodic variables, especially the binary candidates, 
presents difficulties in the confirmation of their nature and further study of their properties. 
We believe, however, that this dataset will be an important resource for variability studies 
of a large and diverse array of objects, serving as a prototype in advance of future large 
synoptic surveys such as Pan-STARRS and LSST. 
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Table 1. Data fields extracted from Stripe 82 object catalogs 



tsObj FITS/CAS Database Field 



Description 



RUN 

RERUN 

CAMCOL 

FIELD 

ID 

RA 

DEC 

OBJC_FLAGS/FLAGS 
TYPE*^ 

PARENT /PARENTID'* 

NCHILD'^ 
STATUS^ 

PSFCOUNTS /PSFM AG 

PSFCOUNTSERR/PSFMAG_ERR 
FIBERCOUNTS /FIBERMAG 
FIBERCOUNTSERR/FIBERMAG_ERR 
MJD 

REDDENING /EXTINCTION 



The SDSS run 

Version of the SDSS FRAMES pipeline (40 or 41 in our dataset) 

CCD camera column (ranges from 1 to 6) 

Field number in SDSS run (13' x 10' per field) 

Non-unique ID number of a detection in a field 

Right ascension J2000 (°) 

Declination J2000 (°) 

Quality flags associated with a detection 

Detection type (3 = galaxy, 6 = star) 

ID of parent object if current detection is a deblended child 

Number of deblended children if current detection is a blend 

Status of the detection (see text for details) 

PSF fitted magnitude in ugriz bands 

Uncertainty in PSF fitted magnitude in ugriz bands 

Magnitude from flux in 3" fiber radius in ugriz bands 

Uncertainty in fiber magnitude in ugriz bands 

MJD of detection of object in ugriz bands 

Extinction in magnitudes for ugriz bands 



'^These fields are used to select suitable objects for extraction and are not saved in the output catalog 
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Table 2. Color selection algorithms for classifying point sources 



Object Type 



Color Selection 



0.98 <{u-g) < 1.30 and -0.05 < RRl < 0.35 and 0.06 < RR2 < 0.55 

and -0.15 < (r - i) < 0.22 and -0.21 < {i - z) < 0.25 

(u-g) < 2.25 and -0.2 < {g - r) < 1.2 and 0.5 < {r ~ i) < 2.0 and 

-19.78(r + 11.13 < {g - r) < 0.95(r ~ i) + 0.5 

-1.0 < (g-r) < -0.2 and -1.0 < (u - g) < 0.7 and 

u- g + 2{g-r) < -0.1 

-0.5 < (5 - r) < 0.75 and 0.6 < (u - 5) < 3.0 and I > 0.135 
2.5 <{u- g) < 3.5 and 0.9 < (5 - r) < 1.3 and s < -0.06 
-1.0 < {u - g) < -0.2 and -0.5 < (5 - r) < 0.2 
0.20 < (.g - r) < 0.48 

-0.7 < PI < -0.25 and 0.4 < {u - g) < 1.4 and -0.5 < (5 - r) < 0.7 
0.48 < (g - r) < 0.55 
0.55 < {g - r) < 0.75 

0.7 <{u- g) < 4.0 and 0.35 < (5 - r) < 0.7 and 

0.15 <{r-i) < 0.6 and I > 0.07 

(.9 - r) > 1.6 and 0.95 < {r - i) < 1.3 

0.666 < 0.875(r - i) + 0.484[(i - z) + 0.00438] < 3.4559 

z < 19.5 and u > 21.0 and g > 22.0 and r > 21.0 and (z - z) > 1.7 



RR Lyrae candidate'^ 

Main-sequence + white-dwarf*^ 

sdO / sdB /white-dwarf'' 

low-metallicity*' 

A/blue horizontal branch'' 

F-turnoff/sdFi' 

Qb 

K-dwarf'^ 
K-giant^ 

sdM'' 

M-dwarf 

brown-dwarf^ 



Using colors from lSesar et al.l (|2007l ). 



''Using colors from lYannv et al.l (2009). 

Color locus computed from mean M-dwarf colors in West et al. ( 20081 ). 
''Using magnitudes that are not dereddened. 



Note. — The various color indices used above are defined below: 
Pl-color = 0.91(?i - g) + 0.415(g - r) - 1.28 
Z-color = -0.436m + 1.129^ - 0.119r - 0.574i + 0.1984 
s-color = -0.249u -I- 0.795g - 0.555r + 0.124 
i?Pl-color ^ {u~g) + 0.67(5 - r) - 1.07 
i?P2-color = 0.45(u ~ 5) - (5 - r) - 0.12 
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Table 3. Distribution of point sources by type and variability (r < 22.0 sample) 



Object type 


Probable variables 


Total objects 


Variable fraction 


sdO /sdB /WD 


4 


235 


0.017 


A/BHB 


180 


1841 


0.097 


F-turnofF/sdF 


507 


26294 


0.019 


F/G 


1333 


42901 


0.031 


G dwarf 


130 


8346 


0.016 


K giant 


141 


14790 


0.009 


AGE 


14 


6001 


0.002 


Low-metallicitv 


348 


13522 


0.026 


K dwarf 


193 


16141 


0.012 


M subdwarf 


4 


332 


0.012 


M dwarf 


253 


107483 


0.002 


brown dwarf 








0.000 


MS+WD 


82 


602 


0.136 


unknown*^ 


1102 


8463 


0.130 


RR Lyrae candidate*' 


125 


1664 


0.075 


SDSS QSO= 


2704 


4196 


0.644 


SDSS white dwarf^ 





419 


0.000 


SDSS hot subdwarf* 





15 


0.000 


SDSS-I variable*^ 


2766 


3972 


0.696 


SDSS standard'' 


562 


165591 


0.003 


All point sources 


6520 


221842 


0.029 



^No color classification possible and no matches to other catalogs. 
^'Using color-selection from lSesar et al.l ( 20071 ). 



'^Cross-matched to objects in lSchneider et al.l (j2007l ). 
^Cross-matched to objects in lEisenstein et al.l (120061 ). 



"^Cross-matched to objects in ISesar et al.l (j2007l ). 
^Cross-matched to objects in llvezic et al.l (j2007l ). 



Note. — Objects may have multiple types assigned, based on their colors. 
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Table 4. Periodic variable candidates classified by light-curve shape 



Periodic variable type Number 

RR Lyrae type RRab 36 
RR Lyrae type RRc 19 
Delta Scuti 16 
Eclipsing Binary 30 



Total 



101 
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Table 5. Recovery rates for artificial sinusoidal variables as a function of r magnitude bin 



Magnitude Bin 


Inserted Objects 


Tentative Variables 


Probable Variables 


Recovery Fraction 


16.0 < r < 17.0 


9 


6 


6 


0.67 


17.0 <r< 18.0 


103 


88 


88 


0.85 


18.0 <r< 19.0 


180 


158 


157 


0.87 


19.0 < r < 20.0 


303 


246 


245 


0.81 


20.0 < r < 21.0 


454 


362 


360 


0.79 


21.0 < r < 22.0 


344 


241 


235 


0.68 



Table 6. Eclipsing and Ellipsoidal Binary Candidates 





Object 


Obs. 


u 


9 


r 




i 


z 


Period [days] 


Comments 


SDSS 


J000845.39+002744.2 


54 


20, 


.51 


19.36 


18, 


.82 


18.60 


18.53 


0.34419 


low-metallicity 


SDSS 


J002719.16+002400.6 


47 


21. 


.04 


20.08 


19, 


.74 


19.61 


19.58 


1.31839 


F/G 


SDSS 


J002851.08+000751.0 


49 


22 


.70 


20.50 


19, 


.07 


18.04 


17.48 


0.59098 


M2 


SDSS 


J003042.11+003420.2 


51 


23, 


.02 


20.57 


19, 


.08 


17.93 


17.32 


0.45705 


M3 


SDSS 


J011155.73-002633.0 


34 


22, 


.37 


20.57 


19, 


.94 


19.71 


19.55 


0.22758 


K dwarf, semi-detached? 


SDSS 


J011156. 52-005221.4 


35 


19, 


.25 


18.29 


18, 


.02 


17.95 


17.99 


0.23002 


F/G, contact? 


SDSS 


J011302. 57+004822.9 


52 


21, 


.20 


19.34 


18, 


.50 


18.14 


17.96 


0.31660 


F/G? 


SDSS 


J011405.02+001138.5 


44 


20, 


.78 


19.88 


19, 


.67 


19.62 


19.64 


0.28550 


F/G, contact? 


SDSS 


J013536.05-011058.7 


44 


20, 


.99 


20.02 


19, 


.95 


19.99 


20.06 


0.39226 


A/BHB, contact? 


SDSS 


J015429.30+005326.7'' 


43 


21, 


.91 


19.50 


18, 


.17 


17.25 


16.77 


2.63902 


Ml 


SDSS 


J015940.01+010328.4 


58 


23, 


.46 


21.47 


20, 


.03 


19.36 


18.97 


0.33779 


MO 


SDSS 


J020540.08-002227.6 


57 


20, 


.39 


19.42 


19, 


.16 


19.10 


19.12 


0.79789 


F/G, deep primary 


SDSS 


J020816.51+003510.0 


48 


19, 


.69 


18.61 


18, 


.77 


18.93 


19.07 


1.61420 


A/BHB 


SDSS 


J021121. 55-003808.3 


40 


21, 


.57 


19.10 


17, 


.77 


17.07 


16.63 


0.31210 


MO 


SDSS 


J021624.34-001817.7 


45 


21, 


.95 


20.80 


20, 


.43 


20.34 


20.28 


0.63678 


F/G 


SDSS 


J022733.94+002615.2 


49 


20, 


.04 


18.83 


18, 


.96 


19.09 


19.21 


0.70626 


A/BHB 


SDSS 


J022858.99-004120.4 


37 


20, 


.61 


19.53 


19, 


.55 


19.63 


19.73 


0.64725 


A/BHB 


SDSS 


J023621.96+011359.1 


57 


19, 


.97 


19.00 


18, 


.74 


18.67 


18.66 


0.35666 


F/G, contact? 


SDSS 


J024109.55-h004813.6 


52 


21, 


.27 


19.29 


18, 


.31 


17.90 


17.64 


0.27645 


deep primary 


SDSS 


J024255.78-001551.5 


62 


21, 


.42 


20.20 


19, 


.68 


19.48 


19.39 


3.19764 


low-metallicity, long period 


SDSS 


J025953.33-004400.3 


53 


21, 


.46 


20.50 


19, 


.36 


18.26 


17.57 


0.14418 


M2, contact? 


SDSS 


J030753. 52+005013.0 


115 


21, 


.84 


19.74 


18, 


.75 


18.34 


18.12 


0.35353 


K dwarf, detached 


SDSS 


J030834.42+005835.2 


68 


22, 


.93 


20.61 


19, 


.49 


19.01 


18.78 


0.27080 


contact? 


SDSS 


J031002.47-000916.2 


54 


19, 


.79 


18.83 


18, 


.44 


18.31 


18.29 


2.19058 


F/G 


SDSS 


J031021.22+001453.9 


48 


19, 


.77 


18.62 


18, 


.11 


17.92 


17.82 


0.26684 


low-metallicity, contact 


SDSS 


J031823.88-010018.4'^ 


34 


22, 


.80 


20.56 


19, 


.10 


17.55 


16.73 


0.40704 


M4 


SDSS 


J032515.05-010239.7 


48 


22, 


.22 


19.86 


18, 


.39 


17.42 


16.87 


0.39451 


M2, shallow eclipses 


SDSS 


J032949. 18-001240.8 


62 


23 


.96 


21.16 


19, 


.78 


19.20 


18.85 


0.39196 


MO 



Table 6 — Continued 



Object 


Obs. 


u 


.9 


r 


i 


z 


Period [days] 


Comments 


SDSS J034256.26-000058.0 


51 


20.29 


19.40 


19.13 


19.06 


19.09 


0.32034 


F/G, semi-detached? 


SDSS J034757.70-001423.5 


53 


20.51 


19.68 


19.32 


19.24 


19.25 


0.27509 


F/G 


SDSS J035138. 50-003924.5 


54 


21.03 


19.07 


18.18 


17.88 


17.70 


0.19892 


K/M?, contact 


SDSS J035300. 50+004836.0 


40 


22.39 


20.62 


19.31 


18.17 


17.61 


0.14855 


M2, ellipsoidal? 



'^This is the confirmed eclipsing binary found by Becker et al. ( 20081) . 
''This is the confirmed eclipsing binary found by Blake et al. 1 20081 ). 



Note. — 



Objects are presented in order of increasing right ascension. 
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Table 7. RRab and RRc RR Lyrae Variables 





Object 


Obs. 


u 


9 




r 




i 




z 




Period [days] 


Type 


0T\00 


J000803, 


.72-010557.9 


49 


20, 


.61 


19, 


A K 

.45 


1 A 

19. 


1 7 

,17 


1 A 

19, 


.11 


1 A 

19. 


,11 


0, 


.64004 




0T\00 


J001031. 


.09+010132.4 


A O 

43 


20, 


.21 


19, 


.06 


Lo. 


07 

,»7 


1 O 
1», 


.81 


18. 


7f? 

,7d 


0, 


.33755 


1 > T 1 , 


bUbb 


J001301, 


.12-003502.6 


27 


19, 


,41 


1 o 

18 


1 

.21 


1 ^7 

17, 


A A 

,94 


1 7 

17, 


.81 


1 7 

17, 


oo 
,82 


0, 


.61226 


KKab 


bUbb 


J001800, 


.72+010424.9 


66 


21, 


.40 


on 
2U, 


A 1 

.41 


on 
2U. 


o^^ 
,2D 


on 
2U, 


oi 
.21 


on 
2U. 


1 n 

,iy 


0, 


.53879 


KKab 


bUbb 


J003119, 


.26+005055.8 


4(5 


20, 


.41 


ly. 


.61 


ly. 


1 o 

,lz 


1 a 

ly. 


no 


1 o 

ly. 


nK 
,Uo 


0, 


.62216 


KKab 


bUbb 


J003209, 


.57-000445.1 


2d 


20, 


.52 


19, 


A A 

.40 


1 A 

ly. 


/I o 

,42 


1 A 

19, 


.48 


1 A 

ly. 


i^O 

,62 


0, 


.20442 


KKc 


bUbb 


J003748, 


.95-000312.0 


A n 
49 


20, 


,74 


19 


TA 

.70 


1 A 

ly. 


A 1 

,41 


1 A 

19, 


oo 
,32 


1 A 

ly. 


O 1 

,31 


0, 


.62441 


KKab 


bUbb 


J004125, 


.63+010859.3 


1 A 

14 


18, 


.65 


17, 


.53 


17. 


OA 

,39 


1 7 

17, 


.41 


1 7 
1 


A O 

,48 


0, 


.25317 


KKc 


bUbb 


J004951, 


.22-000705.8 


2d 


20, 


.85 


1 n 

ly. 


.00 


1 n 

ly. 


7n 

, IV 


1 n 

ly. 


^?7 
.0 ( 


ly. 


,ol 


0, 


.20296 


KKc 


bUbb 


J005253, 


.90-002141.8 


2 / 


20, 


,37 




QA 

.oU 


1 n 

ly. 


no 


lo. 


Arc 

,yo 


1 

lo. 


no 

,y2 


0, 


.58011 


KKab 


bUbb 


J005712, 


.29-000557.7 


2 1 


19, 


,33 


1 o 

18 


O A 

.24 


1 o 

18, 


AA 

,09 


1 o 

18, 


AO 

,08 


1 o 

18, 


1 A 
,10 


0, 


.38487 


TZ> TZ> ... 

KKc 


bUbb 


J012052, 


.63-010411.8 


A O. 

4o 


21, 


,10 


19, 


AO 

.93 


1 A 

ly. 


,08 


1 A 

19, 


,d2 


1 A 

ly. 


,5o 


0, 


.58336 


TD TZ> 1^ 

KKab 


bUbb 


J012333, 


.48+010420.7 


15 


18, 


.89 


17, 


.76 


17. 


CO 

,52 


1 7 

17, 


.46 


1 7 

17. 


,43 


0, 


.28045 


1 > T 1 , 

KKc 


bUbb 


J012543, 


.58+011032.7 


A Q 

4o 


20, 


,51 


1 n 

iy 


A K 
.40 


1 n 

iy. 


A O 

,42 


1 o 

ly. 


,4o 


1 A 

ly. 


C^7 
,0 ( 


0, 


.20111 


KKc 


bUbb 


J012714, 


.38+004403.6 


1 A 

14 


18, 


.79 


1 1. 


^A 

.D4 


1 < . 


,0O 


1 7 
1 /, 


.oU 


1 7 
1 1. 


t^A 

,04 


0, 


.27462 


KKc 


bUbb 


J012813, 


.56+011357.4 


28 


19, 


.12 


17, 


n A 
.94 


17. 


7A 

,70 


1 7 

17, 


.63 


1 7 

17. 


C A 

,59 


0, 


.58792 


KKab 


bUbb 


J013352, 


.87-004502.7 


1 o 

lo 


18, 


.64 


17, 


A T 

.47 


17. 


oo 
,22 


1 7 

17, 


1 A 
.10 


1 7 

17. 


1 A 
,10 


0, 


.60606 


oo ,. 1 , 

KKab 


bUbb 


J014528, 


.54+005133.7 


47 


20, 


,17 


ly 


1 A 

.14 


1 A 

ly. 


,11 


1 A 

19, 


1 7 

,17 


1 A 

ly. 


O A 

,24 


0, 


.22814 


o o ... 
KKc 


bUbb 


J015813, 


.75+010143.5 


6 1 


19, 


.31 


18, 


.1 1 


1 / . 


nn 

,yu 


1 7 
1 ( , 


oo 
.02 


1 7 
1 / . 


7n 


0, 


.70862 


KKab 


bUbb 


J020206, 


.96-003534.2 


O A 

34 


21, 


.97 


on 

20, 


.84 


OA 

20. 


,Dl 


OA 

20, 


c o 
.53 


OA 

20. 


OA 

,39 


0, 


.59791 


oo ,. 1 , 

KKab 


bUbb 


J020245, 


.96-000001.7 


24 


18, 


,61 


1 / 


.51 


1 7 
1 / , 


O/^ 

,2o 


1 7 

1 1. 


1 o 

,18 


1 7 
1 / , 


1 o 

,18 


0, 


.56383 


O O ...K 

KKab 


bUbb 


J020430, 


.98+002005.4 


19 


19, 


,12 


1 / 


A C 

.95 


1 7 
1 / , 


OA 

,89 


1 7 
1 ( , 


o o 
,83 


1 7 
1 / , 


oo 
,82 


0, 


.63040 


KKab 


bUbb 


J020548, 


.43+000144.4 


A A 

40 


19, 


.54 


18, 


A O 

.42 


18. 


OA 

,20 


1 O 

18, 


.14 


1 O 

18. 


,15 


0, 


.55212 


O O «K 

KKab 


bUbb 


J020825, 


.13-003444.4 


6o 


19, 


.68 


Is, 


.0(5 


lo. 


,C>0 


lo. 


Qn 


lo. 


0^5 

,2D 


0, 


.51544 


oo 1 , 

KKab 


SDSS 


J022309 


.05+005734.2 


43 


19, 


,15 


18, 


.02 


17, 


,82 


17, 


,71 


17, 


,70 


0, 


.63008 


RRab 


SDSS 


J022429, 


.81+000725.8 


34 


19, 


.44 


18, 


.38 


18. 


,20 


18, 


.16 


18. 


,16 


0, 


.47630 


RRab 


SDSS 


J022857, 


.91+005359.8 


44 


19, 


.52 


18, 


.47 


18. 


,22 


18, 


.13 


18. 


,13 


0, 


.61506 


RRab 


SDSS 


J023001, 


.99-011146.9 


39 


19, 


.55 


18, 


.43 


18. 


,16 


18, 


.11 


18. 


,11 


0, 


.64910 


RRab 


SDSS 


J023346 


.13-001901.7 


22 


19, 


,07 


17, 


.97 


17, 


,75 


17, 


,69 


17, 


,68 


0, 


.50251 


RRab 


SDSS 


J023635, 


.42+003005.0 


44 


19, 


.65 


18, 


.49 


18. 


,24 


18, 


.17 


18. 


,16 


0, 


.53365 


RRab 


SDSS 


J023932, 


.33+002650.0 


45 


19, 


.53 


18, 


.27 


18. 


,26 


18, 


.31 


18. 


,35 


0, 


.30838 


RRc 


SDSS 


J024216, 


.46+003438.6 


52 


20, 


,30 


19, 


.27 


19, 


,28 


19, 


,35 


19, 


,41 


0, 


.20272 


RRc 


SDSS 


J024647, 


.88+010944.7 


47 


20, 


,36 


19, 


.33 


19, 


,26 


19, 


,32 


19, 


,34 


0, 


.30117 


RRc 


SDSS 


J024701, 


.75-010814.8 


20 


18, 


.72 


17, 


.64 


17, 


,48 


17, 


.49 


17. 


,51 


0, 


.21192 


RRc 


SDSS 


J024727, 


.34-000026.1 


35 


20, 


.26 


19, 


.19 


19. 


,11 


19, 


.21 


19. 


,26 


0, 


.30126 


RRc 


SDSS 


J025658, 


.30-010918.3 


43 


19, 


,09 


18, 


.03 


17, 


,76 


17, 


,70 


17, 


,69 


0, 


.51943 


RRab 


SDSS 


J025717, 


.50+004706.9 


38 


19, 


.79 


18, 


.63 


18. 


,42 


18, 


.36 


18. 


,36 


0, 


.55448 


RRab 
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Table 7— Continued 





Object 


Obs. 


u 


9 




r 




i 




z 




Period [days] 


Type 


SDSS 


J025828, 


.94-003635.1 


27 


19, 


.13 


18, 


.06 


17. 


,91 


17, 


.94 


17. 


,95 


0, 


.35891 


RRc 


SDSS 


J030304, 


.27+003029.6 


53 


19, 


,47 


18, 


.34 


18, 


,11 


18, 


,03 


18, 


,07 


0, 


.64189 


RRab 


SDSS 


J030413, 


.34+005828.1 


54 


19, 


.34 


18, 


.27 


18. 


,00 


17, 


.93 


17. 


,94 


0, 


.61633 


RRab 


SDSS 


J030413, 


.80-011314.3 


19 


18, 


.37 


17, 


.31 


17. 


,19 


17, 


.15 


17. 


,17 


0, 


.36579 


RRc 


SDSS 


J031537, 


.87-005341.9 


30 


19, 


.08 


17, 


.97 


17. 


,74 


17, 


.65 


17. 


,66 


0, 


.62367 


RRab 


SDSS 


J032109 


.31+000705.7 


31 


18, 


,99 


17, 


.86 


17, 


,62 


17, 


,53 


17, 


,51 


0, 


.56928 


RRab 


SDSS 


J032202, 


.45-000151.6 


12 


21, 


.63 


20, 


.49 


20. 


,28 


20, 


.13 


20. 


,09 


0, 


.67613 


RRab 


SDSS 


J032910, 


.97+003614.3 


55 


19, 


.51 


18, 


.32 


18. 


,10 


18, 


.09 


18. 


,11 


0, 


.34381 


RRc 


SDSS 


J032953, 


.69-011413.0 


64 


20, 


,45 


19, 


.46 


19, 


,21 


19, 


,11 


19, 


,11 


0, 


.58634 


RRab 


SDSS 


J033343, 


.92+004457.6 


32 


18, 


,95 


17, 


.84 


17, 


,62 


17, 


,56 


17, 


,59 


0, 


.40888 


RRc 


SDSS 


J033817, 


.21-011112.9 


62 


20, 


.70 


19, 


.69 


19, 


,41 


19, 


,29 


19. 


,32 


0, 


.67096 


RRab 


SDSS 


J034239. 


.98-000009.9 


74 


19, 


.33 


18, 


.25 


17. 


,99 


17, 


.91 


17. 


,93 


0, 


.61181 


RRab 


SDSS 


J034602, 


.52-003342.2 


41 


19, 


,31 


18, 


.20 


18, 


,08 


18, 


,06 


18, 


,13 


0, 


.35783 


RRc 


SDSS 


J034830, 


.07+004042.4 


39 


19, 


.68 


18, 


.57 


18. 


,34 


18, 


.27 


18. 


,24 


0, 


.61444 


RRab 


SDSS 


J034830, 


.94-002320.9 


56 


19, 


.88 


18, 


.79 


18. 


,53 


18, 


.49 


18. 


,50 


0, 


.63686 


RRab 


SDSS 


J034836, 


.34+005334.7 


53 


19, 


.48 


18, 


.42 


18. 


,13 


18, 


.07 


18. 


,08 


0, 


.58038 


RRab 


SDSS 


J034837, 


.58+000804.8 


57 


19, 


,15 


18, 


.13 


17, 


,90 


17, 


,84 


17, 


,87 


0, 


.64671 


RRab 


SDSS 


J035400, 


.50+002110.2 


32 


19, 


.02 


18, 


.01 


17. 


,76 


17, 


.71 


17. 


,73 


0, 


.58990 


RRab 



Note. — 



Objects are presented in order of increasing right ascension. 
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Table 8. High Amplitude Delta Scuti Variables 





Object 


Obs. 


u 


9 




r 




i 




z 




Period [days] 


SDSS 


J001142, 


.02+005946.8 


42 


20, 


,13 


19. 


,17 


19, 


,05 


19. 


,05 


19, 


.09 


0, 


,06837 


SDSS 


J001544, 


.29+003735.1 


84 


20. 


.03 


18. 


,92 


18, 


.92 


18. 


,97 


19, 


.03 


0. 


.04953 


SDSS 


J005312, 


.91-000522.5 


37 


21, 


.14 


20. 


,19 


20, 


.10 


20. 


,13 


20, 


.21 


0, 


.06073 


SDSS 


J013945, 


.70-003325.6 


50 


20, 


.80 


19. 


,73 


19, 


.74 


19. 


,77 


19, 


.82 


0, 


.06480 


SDSS 


J014200, 


.19-001756.6 


49 


20, 


.70 


19. 


,64 


19, 


.58 


19. 


,65 


19, 


.72 


0, 


.06682 


SDSS 


J014225, 


.68+001549.6 


51 


20, 


.27 


19. 


,28 


19, 


.18 


19. 


,18 


19, 


.23 


0, 


.09279 


SDSS 


J014442, 


.67-003741.7 


94 


20, 


.55 


19. 


,57 


19, 


.54 


19. 


,55 


19, 


.54 


0, 


.06245 


SDSS 


J014531, 


,34+011211.9 


56 


20, 


.09 


19. 


,21 


19, 


.10 


19. 


,08 


19, 


.12 


0, 


.08348 


SDSS 


J020841, 


.38-002038.3 


50 


19, 


.60 


18. 


,57 


18, 


.49 


18. 


,54 


18, 


.62 


0, 


.06145 


SDSS 


J021541, 


.11-000750.1 


47 


20, 


.50 


19. 


,55 


19, 


.47 


19. 


,53 


19, 


.57 


0, 


.05419 


SDSS 


J024132, 


.57+002324.1 


37 


20, 


.27 


19. 


,28 


19, 


.18 


19. 


,21 


19, 


.25 


0, 


.06852 


SDSS 


J024751, 


.89-002434.5 


47 


21, 


.13 


20. 


,17 


20, 


.07 


20. 


,10 


20, 


.07 


0, 


.06713 


SDSS 


J030256, 


.28-004150.7 


41 


18, 


.84 


17. 


,86 


17, 


.82 


17. 


,87 


17, 


.93 


0, 


.05813 


SDSS 


J030758, 


.97+000751.7 


58 


20, 


.42 


19. 


,41 


19, 


.41 


19. 


,44 


19, 


.51 


0, 


.06542 


SDSS 


J031956, 


.18-000238.9 


49 


20, 


.00 


19. 


,01 


18, 


.96 


19. 


,00 


19, 


.08 


0, 


.05815 


SDSS 


J034507, 


.01+010040.1 


30 


21, 


.30 


20. 


,16 


20, 


.19 


20. 


,21 


20, 


.13 


0, 


.05733 



Note. — 



Objects are presented in order of increasing right ascension 
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51000 52000 53000 54000 55000 0-0 0.5 1.0 1.5 2.0 2.5 3.0 

MJD log(seporotion[doys]) 



Fig. 1. — Left: Stripe 82 run temporal coverage vs. right ascension. The dense clusters 
at the right are the observation runs associated with the Supernova Survey. Note the large 
gaps between successive years of coverage of the Stripe. Right: Histogram of the separation 
in days between consecutive observations of Stripe 82. The bins are 0.25 dex wide. The 
solid hne represents all observations of the Stripe; the dashed hne represents observations of 
the Stripe before the SN Survey; the dot-dashed line represents observations of the Stripe 
carried out during the SN Survey. The high relative cadence of the SN Survey observations 
is apparent. 
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Fig. 2. — A schematic view of our processing pipeline. The initial extraction, reduction 
of object catalogs from the SN Survey, and ensemble photometry is described in Section [2J 
Variables are detected using the methods presented in Section [3?T1 We find periodic variables 
by running the period search algorithms outlined in Section 13.21 Once these are found, we 
refine their periods, and generate their phased light-curves and ephemerides. Non-periodic 
variables include long term variables (variable on scales of years), quasars, and flare stars. 




Fig. 3. — Left: Histogram of observations per object in the entire dataset spanning the years 
1998 to 2007. The median number of observations is 30. Right: Histogram of observations 
per object for objects with at least 60 observations. 
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Fig. 4. — Top-left: SDSS r median light-curve magnitude vs. SDSS r light-curve standard 
deviation. Top-right: Median SDSS r light-curve standard deviation plotted against median 
SDSS r LC magnitude, binned in 0.5 magnitude bins, for 16.0 < r < 22.0. Bottom-left: 
SDSS r median hght-curve mag vs. differential r hght-curve standard deviation. The scatter 
in the magnitude-cr relation is much reduced. Bottom-right: Median differential r light- 
curve standard deviation vs. median SDSS r LC magnitude, binned in 0.5 magnitude bins, 
for 16.0 < r < 22.0. 'Error-bars' represent the rms scatter for each bin of the respective 
light-curve standard deviation distributions. The lines are fourth-degree polynomial fits. All 
density contours are in half-dex increments ranging from 1.0 to 3.0. 
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Fig. 5. — Left: SDSS r mag light-curve for an M3 dwarf eclipsing binary candidate (plus 
symbols) compared to its four closest neighbors (diamond, square, triangle, and asterisk 
symbols). Note that all five objects appear to have a light-curve dip on the same observation 
date (outlined by the box); this is an example of false variability caused by systematic 
effects from varying photometric conditions. Right: differential r light-curve after ensemble 
photometry carried out on the this field. Note that only the target object (plus symbols) 
appears to be variable: a possible eclipsing binary. The other objects no longer suffer from 
systematic effects that would have resulted in them being erroneously tagged as possible 
variables. 
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Fig. 6. — Examples of PSF r differential magnitude-cx trend plots generated by our ensemble 
photometry pipeline. These particular plots are for the comparison star ensemble fields of 
two objects eventually identified as periodic variables: MB875 (left panel) and MB32047 
(right panel). Objects tagged by the ensemble photometry pipeline as tentative variables are 
shown as star symbols and the two periodic variables in the field are shown as stars inside 
the diamond symbols. 
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Fig. 7. — Stetson variability index histograms for all 365,086 objects with more than 10 
detections in our RA to 4 h light-curve catalog. The dot-dashed lines indicate the standard 
deviations of the distributions, while the dashed line indicates our threshold Stetson index 
value of 0.3. 
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Fig. 8. — Stetson variability index Jg^. vs. SDSS light-curve median r magnitude for all 
365,086 objects with more than 10 detections in our RA to 4 h light-curve catalog. The 
density contours are in dex increments ranging from 1.0 to 3.5, and the dashed line shows 
the minimum threshold value for variability detection, Jgj. = 0.3. 
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Fig. 9. — Left: differential (/-band Dworetsky string length diagram for an eclipsing binary 
candidate (MB6467 in Figure [T8|) . Right: differential ^^-band Stetson string length diagram 
for the same object. The smallest value of the string length in both cases indicates the most 
likely period, in this case; ~0.2668 days. 
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Fig. 10. — Lomb-Scargle periodogram for a typical spectral window function in the range 0.1 
to 100.0 days. Note the strong peaks near 1, 7, and 30 days and their aliases. This window 
function was calculated using the (7-band differential magnitude light-curve from the same 
object as in Figure M The period of this object is marked with the grey dashed line. 
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Fig. 11. — Left: The variability fraction as a function of r median SDSS light-curve magni- 
tude for the r < 22.0 sample of 221,842 objects. Right: The median differential magnitude 
light-curve cr as a function of magnitude for objects marked nonvariable (solid line) and those 
tagged as probable variables (dashed line) for the r < 22.0 sample. 
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Fig. 12. — Left: SDSS g-r/u-g color-color diagram for 221,842 point sources in our final 
light-curve catalog for RA to 4 h. The density contours are in half-dex increments ranging 
from 1.0 to 3.0. Regions A~F indicate positions of various types of objects in color-color 
space: A. white dwarfs, B. low-2; quasars, C. A/BHB stars (including RR Lyrae), D. hi-2; 
quasars, E. main stellar locus, F. faint red objects (such as late M dwarfs and brown dwarfs). 
Right: The same diagram, this time only for 6,520 probable variable point sources identified 
by our pipeline. Note that the QSOs in region B dominate the numbers of probable variables 
detected. 
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Fig. 13. — Phased differential magnitude liglit-curves for a sample of 6 objects with uncertain 
periods, in order of increasing period. The internal designation for an object and its best 
estimated period in days is shown in each panel. Diamond symbols are g light-curve points, 
triangle symbols are r light-curve points, box symbols are i light-curve points and cross 
symbols are z light-curve points. The light-curves are phased around the time of minimum 
light and repeated twice from phase -1 to -f-l to clearly show variability. 
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Fig. 14. — Left: Positions in u — g/g — r color space of 101 periodic variables identified in 
this work. Regions A-F are the same as in Figure [T21 Right: Same as in the left panel, but 
zoomed into region C. 
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Fig. 15. — Top left: Period recovery fraction as a function of the median r magnitude of 
synthetic probable variables. Top right: Period recovery fraction as a function of the number 
of observations. Bottom left: Period recovery fraction as a function of the amplitude. Bottom 
right: Period recovery fraction as a function of the input period. 
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Fig. 16. — A comparison of assigned periods for artificial periodic variables and the period 
recovered by our pipeline. The solid line with unit slope indicates that the recovered and as- 
signed periods are equal. The two dashed curves depict the relation Precovered = ^assigned/ (1± 

Passigned), whilc the dot-dashcd curve depicts the relation Precovered = Passigned/ {Passigned-'^)- 

The solid circles are results of this experim ent for each individual artificial periodic variable. 
Also see Figure 6 in IWatkins et all (120091 ) . 
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Fig. 17. — Period histogram for 30 eclipsing and ellipsoidal binary candidates from this 
dataset, using a bin width of 0.2 days. There is a distinct lack of objects with periods 
greater than 1.0 day. 
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Fig. 18. — Phased differential magnitude ligfit-curves for a sample of 6 eclipsing binary 
candidates identified by our pipeline, in order of increasing period. The internal designation 
for an object and its period in days is shown in each panel. Diamond symbols are g light- 
curve points, triangle symbols are r light-curve points, box symbols are i light-curve points 
and cross symbols are z light-curve points. The light-curves are phased around the time 
of minimum light and repeated twice from phase -1 to to clearly show variability. The 
arrows point out the r-band out-of-eclipse variability observed for MB23368 mentioned in 
the text. 
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Fig. 19. — Left: Fourier parameter R21 calculated for g^-band differential magnitude light- 
curves plotted against the period for RRab (triangles), RRc (boxes), and Delta Scuti variable 
candidates (crosses). The three sinusoidal variable types are easily distinguishable from each 
other. Right: Fourier parameter 02i calculated for g-hand differential magnitude light-curves 
plotted against the period for the same objects in the left panel. 
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Fig. 20. — Phased gri differential magnitude liglit-curves for a sample of 6 RRab RR Lyrae 
candidates identified by our pipeline, in order of increasing period. The internal designation 
for an object and its period in days is shown in each panel. Diamond symbols are g light- 
curve points, triangle symbols are r light-curve points, and box symbols are i light-curve 
points. The light-curves are phased around the time of minimum light and repeated twice 
from phase -1 to +1 to clearly show variability. 
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Fig. 21. — Phased gri differential magnitude light-curves for a sample of 6 RRc RR Lyrae 
candidates identified by our pipeline, in order of increasing period. The internal designation 
for an object and its period in days is shown in each panel. Diamond symbols are g light- 
curve points, triangle symbols are r light-curve points, and box symbols are i light-curve 
points. The light-curves are phased around the time of minimum light and repeated twice 
from phase -1 to +1 to clearly show variability. 
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Fig. 22. — Top left: Period recovery fraction as a function of the median g magnitude 
of synthetic RRab and RRc probable variables. Top right: Period recovery fraction as a 
function of the number of observations. Bottom left: Period recovery fraction as a function 
of the amplitude. Bottom right: Period recovery fraction as a function of the input period. 
The solid hues are for the RRab variables, and the dashed lines are for the RRc variables. 
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Fig. 23. — Distance vs. right ascension distribution of the 55 RR Lyrae in our dataset, 
which is restricted to a right ascension range of 0° to 60°, and a dechnation range of -1.27° 
to +1.27°. 
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Fig. 24. — Phased gri differential magnitude light-curves for a sample of 6 Delta Scuti 
variable candidates identified by our pipeline, in order of increasing period. The internal 
designation for an object and its period in days is shown in each panel. Diamond symbols 
are g light-curve points, triangle symbols are r light-curve points, and box symbols are i light- 
curve points. The light-curves are phased around the time of minimum light and repeated 
twice from phase -1 to to clearly show variability. 
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Fig. 25. — Differential magnitude gri light- curves for a random sample of 16 objects from 
the 2,704 quasars from the SDSS Quasar Catalog (ISchneider et al.ll2007l ) tagged as variables 
by our pipeline. Diamond symbols are g light-curve points, triangle symbols are r light-curve 
points, and box symbols are i light-curve points. The long term variability of these quasars 
is apparent over multi-year baselines. 
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Fig. 26. — Left: Stetson variability index J^i a gainst J„g for 2,704 v ariable QSOs (black 
circles) matched to the SDSS Quasar Catalog (jSchneider et al.l 120071 ) and 1,413 variable 
stellar locus objects (grey circles). The 2,403 new variable QSO candidates identified in 
this work are not included in this plot. Right: Stetson variability index Jg^. against Jug 
for all 6,520 probable variables tagged by our pipeline. The black circles are 2,704 variable 
QSOs matched to the SDSS Quasar Catalog plus 2,403 QSO candidates selected by color 
and variability. The grey circles are the remaining 1,413 variable stellar locus objects. 
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Fig. 27. — g — r /u — g color-color diagram for 11,328 probable nonvariable point sources 
identified by our pipeline. Regions A-F are the same as in Figure [121 The density contours 
are in 0.5 dex increments from 0.5 to 2.0. 



